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The equilibrium constant of the reaction HT+H:O=H.2+HTO involving radioactive 
hydrogen (T) has been determined over the temperature range 289 to 576°K in contact with a 
platinum-charcoal catalyst. The results are summarized in the equations: 


log K =0.292 log 7+336.5/T—1.055 
AF°=+4.83T—1.34T log T—1540 
AH°=0.58T — 1540 
AS° = 1.34 log T—4.25 

AC,° =0.58+0.05 cal./deg./mole 

AH,° = —1540+160 cal./mole 


IT=1.055. 


Comparison of the experimental results with theoretical calculations by W. F. Libby indicates 
the correctness of these latter and confirms the assumption with respect to isotopic species that 
the forces acting within a molecule are unaltered by the change in nuclear mass and that any 
effects on the properties of the molecule are solely due to a different mass moving in the same 


force field. 





HE object of these experiments is the deter- 
mination of the equilibrium constant of the 
reaction HT+HHO=HH+HTO, “‘T”’ being 
“tritium,”’ the hydrogen isotope of mass 3, first 
detected in 1934.1 The equilibrium can be fol- 
lowed by virtue of the radioactivity of tritium? 
which is evidenced by a weak beta-radiation with 
a half-life of 318 years.® 
The interest in this determination lies in the 
effect of isotopic substitution upon the potential 
function of a molecule. The usual assumption 
of this effect is that the forces acting within the 


ae Harteck, and Rutherford, Nature 133, 413 


( 
1939" W. Alvarez and Robert Cornog, Phys. Rev. 56, 613 
(1940) D. O’Neal and M. Goldhaber, Phys. Rev. 58, 574-5 


molecule are unaltered by the change in nuclear 
mass and that any effects upon the properties of 
the molecule are solely due to a different mass 
moving in the same force field. On the basis of 
this assumption, calculations have been made by 
W. F. Libby‘ of the equilibrium constant of the 
reaction investigated in this report. This reaction 
is simple enough for accurate theoretical treat- 
ment and is concerned with a mass difference 
which should provide the best possible criterion 
of the assumptions involved concerning isotopes. 


EXPERIMENTAL METHOD 


The equilibrium constant determined here is 


written (HTO)(HH)/(HT)(HHO)=K. Alter- 


4W. F. Libby, J. Chem. Phys. 11, 101 (1943). 
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nately, this may be put 


(HTO)/(HHO) 
(HT)/(HH) 





which is the form into which substitution was 
made in these experiments. 

The partial pressure of HTO (or of HT) ina 
given container is proportional to the number of 
disintegrations resulting from it per minute. 
In a sample of radioactive water or hydrogen 
which will give a convenient counting rate, how- 
ever, the actual mole fraction of tritium present 
is of the order of 10-" so that the total pressure 
of a sample of radioactive hydrogen or water 
vapor may very safely be used as the partial 
pressure of the inactive species. Therefore, after 
equilibrium has been established, by separating 
the hydrogen, active plus inactive, from the 
water vapor, active plus inactive, and putting 
each in turn into the same counter where the 
counts per minute and the pressure to which the 
sample had been introduced were noted, we 
could write counts per minute per cm Hg for 
water vapor (RHTO/HHO) over counts per 
minute per cm Hg for hydrogen (RHT/HH) to 
obtain the desired constant. 

The first problem was to obtain equilibrium. 
This was done for the runs at 16°C, 20°C and 
25°C by placing the platinum on carbon catalyst 
within a water-jacketed condenser. Constant 
temperature was maintained by running water 
from a thermostat through the water jacket. 
The radioactive water was placed immediately 
on the catalyst with the hydrogen being cycled 
over it and around a closed system until equi- 
librium had been reached. A portion of the sys- 
tem, enclosing a sample of the equilibrated 
hydrogen and water vapor, was now shut off. 
No further change would occur within this 
sample since the exchange of hydrogen isotopes 
between hydrogen and water is not catalyzed by 
glass.5 The two gases were separated by freezing 
out the water vapor with a cold bath, after which 
each could be counted separately. 

For the runs above 25°C, a very active catalyst 
was used. This was kept at constant temperature 


5 Gould, Bleakney, and Taylor, J. Chem. Phys. 2, 362 
(1934). 
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in a vapor bath. The inactive hydrogen was 
bubbled through radioactive water, held at 
constant temperature, and the mixture passed 
over the catalyst, where equilibrium was es- 
tablished. It then flowed on through a sample 
trap where a sample could be closed off and 
analyzed as in the previous apparatus. 

In both these methods it is important that a 
sample of hydrogen and water vapor, in equilib- 
rium, was obtained. Thus, the gaseous equi- 
librium was investigated. To do this it was neces- 
sary to be very careful that none of the water 
vapor was removed by condensation after equi- 
librium had been established. Such an occurrence 
would result in a concentration of tritium within 
the condensed phase, due to the extreme isotopic 
effect between tritium and hydrogen, and a 
constant below the true value would result. 

The counting was carried out in a Geiger- 
Mueller tube counter of one-liter volume. The 
counting mixture was composed of 0.2-cm 
water vapor, 2.5-cm alcohol, from 0.8- to 2.0-cm 
hydrogen and 2.0-cm argon. This mixture was 
always introduced in the same order, the pres- 
sures being measured by a closed-end manometer 
which was utilized as a McLeod gauge to read 
the water vapor pressure when that gas was 
the active material. The tube, so filled, had a 
plateau (region where the count is independent 
of the wall voltage) about 20 volts long with a 
center about 50 volts above the threshold (the 
wall voltage at which counts first can be de- 
tected). In each count, we tried to set the 
voltage so that the tube would be operating at 
the center of the plateau, where reproducibility 
is most easily obtained. 

Initially, it was necessary to check the re- 
producibility of the counting conditions by 
bringing a standard uranium sample to an exact 
position near the counter and noting the number 
of counts per minute it produced in addition to 
the count being determined. This latter was then 
corrected by the factor necessary to bring 
this additional number of counts to a standard 
value. 

This method was used for the first four runs 
at 20°C, which was the first temperature in- 
vestigated. For the rest of the runs, however, 
such a correction was not employed, since it was 
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APPARATUS AND PROCEDURE 


The experiments were performed in two 
different sets of apparatus shown in Fig. 1. 

System A comprises the evacuating devices, 
counter tube, measuring devices and storage 
bulbs. These pieces of equipment were main- 
tained as shown throughout the entire series of 
runs. At the point marked by the inscription 
“to system B” and by the arrow, system A was 
sealed to the apparatus in which equilibrium was 
obtained. 

It was in this latter part of the apparatus that 
two different techniques were involved. The first 
method, depicted as ‘system B-1,’’ was used to 
obtain equilibrium at 16°C, 20°C and 25°C. For 
higher temperatures, the equipment shown in 
“system B-2” was employed. Both of these 
apparatus were joined to system A at the 


“sample trap,” the stopcock labeled ‘1”’ and_ 


divided by the dotted line being common to 
both systems A and B in each case. 

The auxiliary electrical apparatus necessary 
for the use of the counter tube in system A was 
contained in a high speed, Geiger-Mueller count- 
ing unit purchased from Herbach and Rade- 
man, Inc. 


System B-1 


Investigations were first carried out in system 
B-1. The catalyst was contained in the inner 
tube of the water-jacketed condenser (A). The 
radioactive water was introduced into the feed 
system (B) through stopcock 4 by lowering the 
mercury level within the feed system. 

The catalyst, kept in good thermal contact 
with the walls by means of the excess water 
wetting them both, was held at constant tem- 
perature by water from a thermostat. A circu- 
lating pump (C) forced the water through the 
outer jacket of the condenser, returning it 
through the overflow trap (D) where the tem- 
perature was measured by a thermometer. It was 
necessary to keep the temperature of the room 
above that of the catalyst to prevent radio- 
active water from distilling out of that part of the 
system. 

The hydrogen, generated electrolytically, was 
passed over platinized asbestos to remove any 
oxygen impurity. The water was removed from 
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the resulting gas by means of phosphorus pen- 
toxide, after which stopcock 6 gave access to the 
reaction system. Within the latter, the hydrogen 
was cycled for several hours, long enough for the 
isotopic equilibrium to be complete for gas and 
water vapor. 

The cycling was accomplished by means of the 
piston pump (£) and mercury valve system (F) 
shown in the drawing. Each of the two identical 
mercury valves consisted of a bulb connected to 
capillary tubing and containing a small amount 
of mercury. Alternating contraction and expan- 
sion of the volume between the valves was 
produced continually by the piston pump which 
set the mercury columns into oscillation. With 
stopcock 7 closed, this resulted in a transfer 
of gas through the valves from right to left. 
Since the valves formed part of the closed tubular 
system containing the catalyst, a continuous 
transfer of gas around the system was ac- 
complished. 

Stopcock 7 was inserted for use in filling the 
system with hydrogen prior to a run. At such a 
time the sample trap (G) was practically empty 
because of removal of the sample for counting. 
Since the apparatus operated at atmospheric 
pressure, an amount of hydrogen corresponding 
to the sample removed was inserted. The system, 
after initially being filled with hydrogen, was so 
operated that no foreign gases were admitted 
and only the hydrogen removed as a sample in 
the previous run had to be replaced to commence 
a new experiment. 


System B-2 


The equilibrium was studied at temperatures 
higher than that of the room by means of sys- 
tem B-2. 

The hydrogen was generated electrolytically 
at about one liter per hour with 2.4 amperes 
current. The generated gas traversed a spray 
trap (A), and then a liquid-air trap (B) to dry 
it completely. The water collected in the latter 
could be forced out through the capillary tubing 
leading to stopcock 6. 

The gas next traversed solenoid valve (C) 
which separated the active water from the in- 
active water when the system was not in opera- 
tion. Beyond this the hydrogen entered the 
bubbler (£), passing down through the long 
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center tube and bubbling through the glass 
beads, covered with radioactive water, at the 
bottom. The radioactive water was added by 
removing the ground glass tapered plug (F) at 
the bend just above the bubbler. The hydrogen, 
saturated with radioactive water, next passed 
through glass wool in the upper half of the bub- 
bler to remove spray. The entire bubbler, in- 
cluding the spray trap, was immersed in a water 
bath kept at constant temperature in a large 
Dewar flask. 

Traversing solenoid valve (D), the gas was 
introduced into the section heated by the vapor 
bath and containing the catalyst (G). Here it was 
brought to temperature by passing downward 
through six-millimeter tubing in which numer- 
ous indentations had been punched to facilitate 
heat exchange. Rounding the bottom ‘U”’ of 
the tube, the mixture of hydrogen and radio- 
active water vapor passed upward over the 
catalyst where equilibrium was established. 

Passing through stopcock 5, the mixture, now 
at equilibrium, could either be run out at the 
exit tube or sent into the sample trap (#) 
through stopcock 4 which was inserted to pre- 
vent placing stopcock 5 under the high vacuum 
sometimes present in the sample trap. The gas 
left the sample trap via stopcocks 3 and 2, 
stopcock 1 being closed. 


System A 


This part of the apparatus was used as de- 
picted throughout the whole series of measure- 
ments. It contained the counter tube, storage 
bulbs, pressure measuring devices, and vacuum 
equipment. System A was first sealed to system 
B-1 and later to system B-2, stopcock 1 being 
the same in all three figures. 

The counter tube itself (A) was constructed 
from two-inch o.d. brass tubing one-sixteenth 
inch thick, and had glass ends. The glass and 
metal were sealed together with picene wax to 
give a vacuum-tight container with a volume of 
about one liter. 

Supported very close to but not touching the 
counter tube wall were the sensitivity holders 
(B). These held two sample tubes of uranium 
acetate, kept firmly in place in the tube by a 
layer of paraffin. These tubes were always 


placed each in the same holder and in exactly 
the same position as indicated by an arrow on the 
cork. Thus, théy should always give the same 
number of counts per minute in addition to those 
being measured inside the counter if other condi- 
tions have been reproduced. Their purpose was 
to check this latter fact. 

Gas was admitted to the counter tube through 
stopcock 2. Stopcock 1 led to the sample trap 
while stopcock 3 led to the manometer (C). 

This latter was a closed tube type which 
could be pulled down by opening stopcocks 12, 
13, and 15 to the floor pump. The manometer was 
used in filling the counter to read the pressures 
of the gases directly or by the McLeod gauge 
method. This last was necessary with the radio- 
active water which was expanded to a pressure 
of about two millimeters on being introduced into 
the counter. 

Stopcock 4 led to the storage bulbs. Sealed to 
stopcock 5 was a test tube (£) of liquid absolute 
alcohol, to stopcock 6 a small bulb (F) of liquid 
inactive water, to stopcock 7 a five-liter bulb 
(G) for inactive hydrogen, to stopcock 8 a liter 
bulb (#) for calibration purposes and to store 
radioactive hydrogen, and to stopcock 9 a five- 
liter bulb (J) of argon. After separation of a 
sample, the radioactive water vapor was stored 
in the sample trap which is not shown as part 
of system A. 

Stopcock 10 led to the main vacuum line. A 
low pressure McLeod gauge (J) tapping this 
line was used to test the efficiency with which 
the apparatus was being evacuated. The upper 
end of the main vacuum line was connected to a 
trap which led into the mercury vapor pump 
(K) backed by a floor pump. 

The vacuum system was in continuous opera- 
tion, since the counter tube had to be evacuating 
at all times when not in use and for 25 minutes 
between each individual counting period. Un- 
less the tube were treated to this exhaustive 
evacuation, the sensitivity would not, in any 
sense, be reproducible but would depend upon 
the recent history of the apparatus. The pe- 
culiarities in evidence were complex and have 
not been explained, but the above treatment 
enabled these measurements to be made with a 
sensitivity more constant than it was possible 
to measure. 
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CALCULATION OF RESULTS 


The data which entered immediately into the 
equilibrium constant were the pressures of HTO 
and HT samples, the counting rate with these 
samples and the background counting rate. 
In the first four runs the counting rates were 
corrected to a standard sensitivity value of 
one count per minute. 

The constant recorded for each temperature 
is the average of the successful runs. The con- 
stant at 25°C, calculated in this manner, is 
6.25+0.05. The error was obtained from the 
actual deviations of the individual results from 
the mean. (See Table I.) 


INTERPRETATION OF RESULTS 


The results of these measurements and their 
conversion into the form necessary for a plot of 
log K against 1/7 are given in Table II. The 


TABLE I. Data and results. 

















HT Obs. HTO 
Temp. Back- Obs.H Press.” —— HTO Press.s -——— 
Page °C ground Count H cm Water Water cm K 
Vol. 2 
14 20.2 0.560 1.873 146 0.900 2.002 0.251 5.74 6.38 
17 20.2 0.583 1.880 1.46 0.888 1.752 0.218 5.36 6.04 
20 20.2 0.560 1915 141 0.961 1.715 0.192 6.01 6.25 
23 20.2 0.567 1.903 146 0.914 1.658 0.199 5.49 6.01 
94 20.2 622 1713 1.45 752 1602 0.179 5475 7.28 
97 20.2 626 1846 =—-:1.48 824 1584 0.181 5300 6.44 
100 20.2 610 1795 =—-1.44 823 1556 80.176 5375 ~~ 6.54 
103 20.2 633 1830 =1.47 814 1572 0.176 §=65340~—6..55 
106 20.2 612 1785 1.48 792 1546 §=6©0.168 5560 7.01 
109 20.2 622 1889 . 1.50 844 1536) «=(0.176 3=— 5200 -~—Ss«6..16 
Av. Ky.2=6.47+0.12 
29 16.0 580 1900 §=1.48 891 1x17 0.140 §=5980_~—s 6.71 
32 16.0 588 1900 §=1.56 841 1323 0.185 5440 6.46 
35 16.0 588 1917 1.51 880 1236 «60.140 §=(5840 (6.63 
38 16.0 600 1932 1.57 848 1404 0.140 5745 6.77 
46 16.0 603 1764 = 1.82 879 1460 0.142 6035 6.86 
49 16.0 605 1725 1.31 855 1440 0.142 5880 6.88 
54 16.0 591 1663 = 1.27 843 1428 0.148 5655 6.71 
58 16.0 596 1690 1.30 841 1367 0.137 5630 6.69 
61 16.0 583 1688 1.32 837 1383 0.140 5710 6.82 
64 16.0 590 1702. 11.33 836 1396 §60.143 §=5640 ~=—-6..74 
67 16.0 585 1640 =1.28 825 1386 §=60.140 §=5720 «6.93 
Av. Ki1s.0o=6.75+0.04 
70 25.0 612 2338 1.99 866 1906 0.238 5440 6.28 
73 25.0 616 1824 1.43 845 1894 0.244 5240 6.20 
76 25.0 600 1860 =1.43 880 1968 0.256 5340 6.07 
79 25.0 604 1862 1.45 867 1875 0.233 5450 6.28 
82 25.0 620 1955 =1.47 907 2199 0.282 5600 6.17 
85 25.0 603 1896 1.43 884 2082 0.262 5640 6.38 
fatal 25.0 623 1888 1.44 878 2051 0.256 5580 6.35 


Av. Ks.0=6.25+0.05 





124 56.3 705 «1889 1.01 1172 1698 0.169 5870 5.01 
127 56.0 625 1860 0.96 1287 1659 0.166 6230 4.84 
133 56.0 665 100 1429 1739 0.151 7140 4.96 
142 56.2 645 2251 1.01 1590 2168 0.196 7770 4.89 
145 56.3 669 2036 0.94 1454 1778 0.151 7340 5.04 
148 = 56.1 675 =621381) = 1.00) 1456 )Ss 1802s 0.151 §=— 7460 = 5.12 
151 55.8 634 = 2153 (0.99) :1533) «1893 0.158 = 7960 = 5.19 
154 55.8 665 2282 1.01 1600 1886 0.146 8360 5.22 
156 55.5 720 0«=—-2333) 1.02) :11582)s«1991 = 0.156 = 8150S 5.15 


Av. Koe.2=5.05+0.05 
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TABLE I.—Continued. 














HT Obs. HTO 
Temp. Back- Obs.H Press) -——- HTO Press. 

Page °C ground Count H em Water Water cm K 
160 79.8 645 2697 0.99 2073 2500 «0.213 «698700 = 4.20 
162 79.8 631 2474 390.86 2144 2309 «490.186 9020 4,21 
165 79.6 631 2400 =—0..81 2184 2335 0.181 9410 4.52 
167 79.5 646 2537 = 0.94 2012 2249 0.173 9260 4.61 
169 79.5 663 2524 30.88 2116 2229 «#20.171 9150 4,32 
171 79.2 642 2532 = 0.89 2124 2603 0.208 9425 4.44 
174 79.6 643 2622 0.94 2105 2445 0.198 9090 4.32 
175 79.7 643 2452 (0.88 2055 2506 0.208 8955 4.36 


Av. K7.0=4.37+0.05 








177) 111.5 676 2585 0.90 2120 «2198 (0.196 967760) = 3.66 
178 111 673 =. 2460 = 0.89 = 2007) 2201) = (0.196) 7795S 3.88 
180 111 673 =. 2367 ~—s «0.81 2091 2069 0.183 7625 = 3.65 
182 111 671 2334 0.80 2080 2325 0.198 8350 4.01 
183111 671 2537 0.83 2248 2190 0.186 8160 3.63 
185 110.0 660 2217 0.80 2217 2255 0.201 7930 3.82 
2 111.2 657 2357) «0.82 2073)=S 2306 = 0.208 )=— 7925 = 3.82 
4 111.3 657 2365 0.82 2084 2254 0.211 7560 3.63 
6 111.2 671 2484 0.83 2185 )=S 2618 )=-0.238 = 8180) 3.74 


Av. Kin.2=3.76+0.05 








7 158.7 656 2389 =—.0.82 2115 1973 =0.191 6900 = 3.26 

9 158.1 656 2332 ~=—-0..80 2095 1980 §$0.196 6750 3.22 
10 =:158.2 656 2358 0.82 2076 1994 0.198 6750 3.25 
11 158.2 685 2494 0.81 2232 2020 = 0.201 6640 = «2.98 
13 158.7 656 2384 390.78 2215 2060 §80.226 ©6210) 2.81 
15 158.6 656 2542 =0..80 2358 2178 = 0.211 7210) = -3.06 
17 158.4 686 2466 0.80 2225 2121 0.204 7030 3.16 
18 158.7 686 2532 0.81 2280 2186 0.208 7200 3.06 

Av. Kiss.4=3.10+0.06 

19 216.9 680 2376 = 0.77 2200 1903 = 0.204 6000 2.72 
21. =217.3 680 2483 = 0.80 2255 1858 0.204 5770 2.56 
22 217.3 680 2450 3=60.81 2185 1836 60.198 5990 2.74 
232s. 217.3 680 2483 =: 0.80 2255 1756 «0.186 5785 2.56 
24 «217.6 664 2653 «(0.91 2185 1654 0.176 5625 2.58 
26 0 8=. 216.5 664 2394 30.77 2247 1815 0.193 5960 2.65 
27 = 217.0 664 2495 0.80 2288 1820 0.181 6390 2.79 
28 = =217.2 664 2643 = (0.81 2442 1888 0.201 6090 2.50 




















30 3023 660 2545 0.81 2327 1714 0.201 5250 2.26 
31 302.3 660 2553 0.81 2338 1720 0.209 5070 2.17 
32 303.0 660 2499 0.80 2287 1638 0.188 5200 2.27 
33 302.9 660 2600 0.83 2336 . 1620 0.193 4975 2.13 
35 303.6 660 2690 0.86 2360 1611 0.186 5110 2.16 
36 6303.6 «= «660 Ss «2530S 0.79 Ss 2366 = :1578 0.186 §=— 4935 2.09 
37 302.8 660 2698 O81 2515 1642 0.186 5280 2.10 
38 302.1 660. 2605 0.81 2400 1627 0.186 5200 2.17 
Av. Kyo2.9=2.17+0.02 
TABLE II. 

FC 7"s. 1/T K log K + 

16.0 289.1 0.003459 6.75+0.04 0.8293 0.0025 
20.2 293.3 3409 64740.12 0.8109 0.0080 
25.0 298.1 3354 6.2540.05 0.7959 0.0034 
56.2 329.3 3036 5.05+0.05 0.7033 0.0043 
79.6 352.7 2834 4.3740.05 0.6405 0.0049 
111.2 384.3 2602 3.76+0.04 0.5752 0.0051 
158.4 413.5 2317 =3.10+0.06 0.4914 0.0084 
217.1 490.2 2039 2.64+0.04 0.4216 0.0065 
2.9 1736 2.17+0.02 0.0040 


576.0 


0.3365 








plot itself appears in Fig. 2. Such a plot is 
usually a straight line but, in this case, due to 
the wide temperature range covered by the ex- 
periments, a slight curvature is present. 

A plot of the slope against temperature is 
utilized here to obtain an estimate of AC,. The 
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data read from Fig. 2 and the slopes obtained 
from them are found in Table III. The slopes are 
plotted against temperature as the lower series 
of points in the bottom part of Fig. 3. The upper 
series of points in this section of Fig. 3 represent 
the deviations of the lower series from the 
straight line drawn through them. From the two 
straight lines that are drawn we find the best 
equation for the data to be slope = —0.126T 
+336=(d log K)/d(1/T). -Multiplying through 
by 2.303R, we obtain 


d log K 
03R———— = — 0.587'+ 1540. (1) 

d(1/T) 
From Eq. (1) we see that AC, is 0.58 cal. per 
mole and A/Z,° is — 1540 cal. per mole. Therefore 
AH®° =0.587 — 1540. (2) 


Integration of Eq. (1) gives us 


—- 336.5 
log K =0.292 log alae iq, (3) 


EQUILIBRIUM 
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Substituting into Eq. (3) the experimental 
values for K with their corresponding tempera- 
tures, we obtain the values for “J’’ shown in 
Table IV. These values show no drift over the 
300° temperature range. We now have the 
complete equation 


~ 


336.5 
log K =0.292 log 7+———— 1.055. (4) 
T 


From the values of AC, and “J’’ we obtain 
for the entropy 


AS° = 1.34 log T7—4.25. (5) 


Combining this with our previous expression 
for AH® we have: 


AF°= —4.83T —1.34T log T—1540. (6) 


To test our equations we use them to calcu- 
late K and compare the value obtained, Keate, 
with that of our experiments, K.x,, and with 
that, Kineo, calculated from pure theory by W. 
F. Libby. This comparison is contained in 
Table V. In each case our calculated value of K 
lies within the root mean square error of the 























TABLE III. 
1/T °K log’K 
log ’K log K 1/T’ 1/T Av. T°K —-—logK Slope 
0.4356 0.3293 0.0021 0.0017 0.0019 526 0.1063 267 
0.4630 0.3546 22 18 0.0020 500 0.1084 271 
0.4910 0.3816 23 19 21 476 0.1094 276 
0.5200 0.4088 24 0.0020 22 454 0.1112 278 
0.5486 0.4356 25 21 23 435 0.1130 282 
0.5764 0.4630 26 22 24 417 0.1134 284 
0.6050 0.4910 27 23 25 400 0.1140 285 
0.6344 0.5200 28 24 26 385 0.1144 286 
0.6630 0.5486 29 25 27 370 0.1144 286 
0.6928 0.5764 0.0030 26 28 357 0.1164 291 
0.7224 0.6050 31 27 29 345 0.1174 293 
0.7520 0.6344 32 28 0.0030 333 0.1176 294 
0.7810 0.6630 33 29 31 323 0.1180 295 
0.8105 0.6928 34 0.0030 32 313 0.1177 294 
TABLE IV. 
—0.292 
T°K log Kos log T logT —336.5/T I 
289.1 0.8293 2.461 0.719 1.164 — 1.054 
293.3 0.8109 2.467 0.720 1.148 — 1.057 
298.1 0.7959 2.474 0.722 1.129 — 1.055 
329.3 0.7033 2.517 0.734 1.022 — 1.053 
352.7 0.6405 2.547 0.744 0.954 — 1.057 
384.3 0.5752 2.585 0.755 0.876 — 1.056 
431.5 0.4914 2.635 0.769 0.780 — 1.058 
490.2 0.4216 2.690 0.785 0.687 — 1.050 
576.0 0.3365 2.760 0.806 0.585 — 1.054 
Av. I= —1.055+0.002 
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experimental value, with no drift being shown 
in these slight deviations. 

The relationship between the experimental 
constants and W. F. Libby’s theoretical calcula- 
tions is demonstrated in Fig. 4. In it the lower 
curve represents the theory and the upper curve 
these experiments. The radius of the circles 
represents, in the case of the experiments, the 
root mean square error of the point while the 
circles around the theoretical points are an 
estimate by W. F. Libby of their probable 
errors. It is interesting to note that in a direction 
perpendicular to the curves, the ‘‘smearing”’ 
produced by these errors just overlaps. 

The small difference between the two curves 
corresponds at 16°C to a AF® difference of 24 cal. 
per mole out of about 1100 cal. This is much too 











TABLE V. 
£°C log Kealc Keale Kexp Ktheo 
16.0 0.828 6.73 6.75+0.04 6.47 
20.2 0.813 6.50 6.47+0.12 6.24 
25.0 0.796 6.25 6.25+0.05 6.01 
56.2 0.701 5.02 5.05+0.05 4.84 
79.6 0.643 4.40 4.37+0.05 4.23 
111.2 0.576 3.77 3.76+0.04 3.64 
158.4 0.494 3.12 3.104:0.06 3.03 
217.1 0.417 2.61 2.64+0.04 2.54 
302.9 0.336 2.47 2.17+0.02 2.08 








high to be due to error in the fundamental 
frequencies of the water molecule. The dis- 
crepancy, therefore, must be due either to some 
unaccounted source of systematic error in the 
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experiments or to some inadequacy in the as- 
sumptions used in the theoretical calculations. 
One of the possibilities in the latter case, and the 
initial interest behind W. F. Libby’s calcula- 
tions, is that the potential function for a molecule 
may, in fact, be changed by isotopic substitution. 

The agreement is actually so close, however, 
as to constitute quite a satisfactory check over 
the entire temperature range. It can be said 
with confidence that isotopic substitution by 
elements other than hydrogen will produce no 
effect upon the potential function of a molecule. 
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A study was made of the exchange of oxygen and carbon 
in the reactions: 


COz (dissolved) +H:O @H2COs, | 
COz (dissolved) +OH-~@HCO;-. 


It was found that the rate of the first could be increased 
four or five times and the rate of the second up to 2000 
times by heterogeneous catalysis. Investigation of the use 
of these catalyzed reactions for separation of the oxygen 
and carbon isotopes was carried out both by the use of a 
static bomb technique and by the operation of fractionating 
columns. The effects of different operating conditions were 


indicated, and the necessity for turbulent flow in the gas 
phase was established. Pressures from one through fifty 
atmospheres were tried, as the increase of CO: pressure 
increases the concentration of COs» (dissolved) and thus 
speeds up the exchange. 

A comparison of the above method with the HCN —CN- 
method of separation of carbon isotopes shows the cyanide 
exchange to be somewhat faster and to require less time 
to come to production, somewhat doubtful advantages 
when they must be coupled with the polymerization and 
poisonous properties of the hydrogen cyanide, as well as 
with its greater cost. 





INTRODUCTION 


NE comparatively new but, nevertheless, 
powerful tool in the study of carbon and 
oxygen is the use of isotopic tracers, which also 
show promise of value in certain types of quan- 
titative analysis. Chemical methods of obtaining 
small quantities of material with different iso- 
topic concentrations than normal have been used 
for some time. 

The two-phase system, sodium cyanide solu- 
tion—hydrogen cyanide gas was successfully 
used! to concentrate C™, but possessed the dis- 
advantages of polymerization of the hydrogen 
cyanide in the presence of sodium cyanide and 
of the highly poisonous character of hydrogen 
cyanide. The method used here was the operation 
of a packed fractionating column by passing 
hydrogen cyanide up the column countercurrent 
to the sodium cyanide solution, and supplying 
reflux at the top. Similar systems using carbon 
dioxide in exchange with bicarbonate ion in 
solution were tried with many types of homo- 
geneous catalysts but without success,? as the 
exchange was very slow and not suitable for 


* Submitted in partial fulfillment of the requirements 
for the degree of Doctor of Philosophy, in the Faculty of 
Pure Science, Columbia Universitv. 

t Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 
— Stewart, and Urey, J. Chem. Phys. 8, 523 

). 
og) Mills and H.C. Urey, J. Am. Chem. Soc. 62, 1019 
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these counter-current methods. We decided that 
investigation of the possibilities of heterogeneous 
catalysis for that exchange might yield an 
exchange rate high enough to be suitable. The 
reactions between carbon dioxide and water 
result in a fractionation factor favorable for an 
isotopic concentration of both C™ and O'%, 
respectively : 


C8O.+xH2C"O3+ (1—x)HC"O3- 
= C"0.+xH2C"0;+4+ (1—x)HC®O3- 
CO.!*+ 2H,0!8= CO,'8+ 2H20"*. 


Since in the cyanide system and other similar 
systems the exchange was much faster than here, 
it follows that the reaction is not diffusion con- 
trolled, but that the slow process is in the 
exchange indicated by (1) and (2): 


k 
COs (dissolved) ++-OH-2HCO;-, (1) 


-1 


k: 
CO, (dissolved) +H.0 = H2CO; (2) 


—2 


Faurholt® in 1924, Brinkman, Margaria, and 
Roughton‘ in 1933, Stadie and O’Brien’ in 1935, 


3C. Faurholt, J. Chim. Phys. 21, 400-455 (1924) (rate 
const. must be multiplied by 2.3). 

4 Brinkman, Margaria, and Roughton, Phil. Trans. Roy. 
Soc. A232, 65 (1933). 

5 W. C. Stadie and H. O’Brien, J. Biol. Chem. 103, 521 
(1933); 112, 723 (1936). 
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Roughton and Booth® in 1938, and Mills and 
Urey’ in 1940 have all determined rate constants 
for some of the above reactions by various 
methods. It is evident that the problem to be 
solved to make possible a practical method of 
using the above exchanges for isotopic concen- 
tration is to increase the rate of the slow process 
while keeping a conservatively low temperature, 
as the fractionation factor decreases with rise in 
temperature.’ 

Cohen® set up equations corresponding to the 
following for isotopic exchange in packed 





columns: 
ON AN 
H—=L—-—kcC{N(1—n)—an(i-N)}, (3) 
Ot 02 


on On 
h—= —l—+kcC{| N(i—n)—an(1i-N)}, (4) 
Ot 02 


where & is a constant; C is the concentration in 
moles per unit volume of carbon in the liquid; 
N is the mole fraction of heavy carbon in the 
liquid ; 77 is the moles of holdup per unit length 
of combined light and heavy carbon in the 
liquid; L is the flow of light and heavy carbon 
in the liquid in moles per unit time; and where 
lower case symbols designate the corresponding 
quantities in the gas phase, z is the distance along 
the column measured from the reflux end, and at 
equilibrium: 


_N(i~n) 
e n(1—N) 


— 
wn 
~~ 


Inspection of (4) and (5) indicates at once 
that for any given column the highest frac- 
tionation is obtained with the highest value of 
kcC. However, that factor is an exchange rate 
constant per unit length of column, and as such 
varies from column to column with the cross- 
sectional area, whereas we wish to find the 
exchange rate per unit volume of packing, irre- 
spective of the size of the particular column used. 
That factor then turns out to be kcC/A, where 


6 F. J. W. Roughton and V. H. Booth, Chem. and Ind. 
56, 1053 (1937). 
' 43 C. Urey and L. J. Greiff, J. Am. Chem. Soc. 57, 321 
(1935). 
8K. Cohen, J. Chem. Phys. 8, 588-597 (1940). 
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A is the cross-sectional area, and is the quantity 
we wish to maximize. 

Introducing the boundary condition that at 
z=0, N=n, and solving (4) and (5), we find that 
for total reflex and no production at z=0; that 
is, /=L, 


kcCz2 2 
—— =——[tanh™ (1—2N.,) 


L a-1 
—tanh (1—2Npo) ]. (6) 


If L=l—P, where P is the number of moles 
per unit time drawn off at the reflux end, then 
LN-In=PNo, andif p=P/L (7) 

kcCz 2 2AN.+B 


- on anh canes 
L(i—p) (B°-—4AD)? (B?—4AD)? 





2AN)+B 
eee (8) 
(B?—4AD)} 


—tanh-! 





where 
A =a- 1, B= 1 —p+pNo—apNo-—a, D=apNyo. 


In these equations, kcCz/L corresponds to the 
number of theoretical plates in the column. 

Corresponding equations are applicable to 
oxygen fractionation. 

Another aspect to be considered is the time 
required for any given column to approach equi- 
librium, as it would be easy to fashion a system 
in which this would take a prohibitive time. 

Cohen finds that when 1/h>L/H and h+Ha 
>H(1—a) 


H 
T= -, 
(1-—a)*LvF(c) 





(9) 


where v is the number of plates/unit length, 7 
being defined as the time required for the 
principal transient term in his expansion to 
diminish to 1/eth of its initial value. F(c) is a 
function of the ratio of the abundance ratios at 
the two ends of the column, decreasing very 
rapidly with increase of that ratio. 

In the case of high pressures in the system, (9) 
would not hold as well, but by multiplication by 
(h+H)/H, one would expect fair agreement, as 
the amount of concentrate to be stripped is 
greater by that factor than in the case where /: 
is negligible. 
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ISOTOPIC EXZCHANGE OF COs, AsCOs, HEOs- 
TABLE I. 
COBO! /— 
(= : ~COx6 Bomb kcC/A 
Exp. Time of . Moles Moles volume moles /cc 

no. reaction Sample Equilibrium T} COs HO in ce min. 
9 5 min. 1.11 1.27 5.62 min. 0.0060 0.0833 Suz 4.06 Xx 1074 
9 15 min. t22 1.27 5.61 min. .0060 0833 3.2 4.05 10~* 
25 20 sec. 1.07 1.44 80 sec. 0.0031 0.0566 5.0 5.83 10-4 
25 30 sec. 1.09 1.44 90 sec. .0031 .0556 5.0 5.18 1074 
EXPERIMENTAL fractionation column to investigate each type of 


Determination of a 


In order to interpret the results of column ex- 
periments properly it is necessary to know the 
value of a. Urey and Greiff? in 1935 calculated 
the constant 


(C0. ][C#0,7] 


of = . ; : 
(C80, J[C"037 ] 


at 25°C to be 1.012. This should be an approxima- 
tion to the value of 


[C02 ]THC"0,-] 
[C#0,]THC®0,;-] 





Urey, Aten, and Keston® found an experimental 
value of about 1.013 or 1.014 for a using the 
equilibrium between KHCO; solution’ and COs. 
In an independent determination we bubbled 
normal CO; for 24 hours through 2 14 NH,HCOs. 
Samples were taken and analyzed on a Nier type 
mass spectrometer until they checked to within 
0.1 percent. A correction was made for the O'7 
abundance, and a was calculated. 





Run Sample a 
A 1 1.0070 
A 2 , 1.0120 
B 1 1.0159 
B 2 1.0132 
Mean 1.0120+0.0025 


The fairly large spread of these values was 
due to spectrometer errors. 

Throughout the calculations in this paper, 
the value of a used was 1.012. 


Static Experiments 


Due to the inconvenience and large time con- 
sumption necessary to set up and operate a 


* Urey, Aten, and Keston, J. Chem. Phys. 4, 622 (1936). 








packing under each set of experimental condi- 
tions, as well as because of the large amount of 
packing necessary for such experiments, we 
decided to investigate first the rates in small 
reaction vessels. Such a vessel would be filled 
with the packing, moistened with the liquid 
phase, filled with the gas phase and allowed to 
stand for a limited length of time. Then analysis 
of the resulting contents would indicate the rate 
of reaction. In such a system, 


j=L=0 hn+HN=hn,t+lIN,. 
If m is small, (4) and (5) yield: 


and 


h H 
0.693—x— 
kcC A A 


=, (10) 


A h+aH 
& 
A 
where 7; 


, is the time required for the system to 
come half-way to equilibrium, and 4/A and J//A 
measure the holdup per unit volume in the gas 
and liquid phases, respectively. 

Copper and steel bombs were loaded with 
packing to be investigated for possible catalytic 
effect. In most cases about 20 percent of the 
total bomb space was filled with water or bicar- 
bonate solution with a lower than normal O'8 
content. This was enough totally to wet the 
packing without apparent surplus. In some cases 
an aqucous solution of bicarbonate with a greater 
than normal C™ content was used. In these cases 
a larger proportion of solution, up to 50 percent 
of the total space, was used so that there would 
be sufficient excess C™ to give a measurable 
enrichment of the added COs. Then the packing 
was at about the limit of liquid holdup without 
flooding. In the latter case the expected exchange 














































REID AND H.C. UREY 
TABLE IT. 
keC/A X105 
Pressure moles/cc min. 
Exp. Temp. Oz 
no. Packing Solution ae. Ib./in.? O- exchange © bs exchange 
1 BaCO; H,0, 0.111% H:O' 23 265 2.77 
2 20% BaCO; on Al,O; H.O, 0.111% H,O'8 23 265 4.04 
3 20% BaCO; on Al;O; H,O, 0.111% H.O'8 23 765 14.8 
4 Diatomaceous earth H;0, 0.111% H:0"8 26 765 14.1 
5 Silica gel H.O, 0.111% H,0'8 24 765 35.0 
6 Calcined silica gel H.O, 0.111% H.O'8 25 765 59 
7 20% PbCO; on Al.O; H2O, 0.111% H20'8 23 85 7.5 
8 20% PbCO; on Al.O; H2O, 0.111% H.O%8 23 265 18.2 
9 20% PbCO; on Al,O; H,0O, 0.111% H20'8 23 765 40, 
10 20% PbCO; on Al,O; H.O, 0.111% H.0'8 0 765 0-2 
11 20% PbCO; on Al.O; H.O, 0.111% H.O' 64 765 80 
12 Al.O3, 15.2% H2O H.O, 0.111% H.0'8 24 265 38 
13 Al.O3, 15.2% H:O H.O, 0.111% H20'8 28 765 140 
14 Al.O3, 15.2% H2O 2.1 M KHCO;, 4.7% KHC#O; 29 765 260 
15 Morell Char. no. 7 H.0, 0.111% H,0'8 13.2 765 127 
16 Morell Char. no. 7 H.O, 0.111% H.0"8 23.5 765 351 
17 Morell Char. no. 7 H.O, 0.111% HO 48 765 351 
18 Morell Char. no. 7* H,O, 0.111% HO 25 765 77 
19 Morell Char. no. 7 2.1 M KHCO;, 4.7% KHC#O; 25 765 210 
20 Morell Char. no. 7 2.1 M KHCO;, 4.7% KHC#O; 25 265 110 
21 Morell Char. no. 7 2.1 M KHCOs;, 3.1% KHC#O; 23.4 72 30 
22 Fiber A: Gassing Fiber 
Glass no. 9930 2.1 M KHCO;, 4.7% KHC*0O; 27 $3 262 
23 Fiber A: Corning Fiber 
Glass no. 9930 2.1 M KHCO;, 4.7% KHC0O; 27 265 450 
24 Fiber A: Corning Fiber 
Glass no. 9930 2.1 M KHCO;, 4.7% KHC#0O; 26 765 630 
25 Fiber B: Corning Fiber 
Glass no. 790 H,20, 0.111% H.0" 23 265 55.0 
26 Fiber B: Corning Fiber 
Glass no. 790 0.50 M KHCO;, 
0.113% [H,08+ KHCO,0'* ] 24 265 123 
27 Fiber B: Corning Fiber 
Glass no. 790 2.1 M KHCO;, 
0.119%[H.0"+ KHCO,0"*] 24 265 280 








* After ageing 41 days. 


TABLE III.* Calculation of constants from kcC/A = (ki [OH~]+2[H20])[CO2 si 




















kcC/A ki [OH-] 

Exp (moles -liter=! +k2[HO] ki (mole - k2[H:0] 
no. Packing *min.~!) [CO] {[OH-] Ww (min.~) min.~!) (min.~) 
22 Fiber A 2.62 0.152 1.311077 2.0 34.8 

2.08 x 108 6.0 
23 Fiber A 4.50 0.63 3.92 x 10-8 2.0 14.3 

2.02 x 108 6.4 
24 Fiber A 6.30 1.18 2.141078 2.0 10.7 

2.05 K 108 6.2 
Average Fiber A 
25 Fiber B 0.55 0.63 1.5x10-" 5.0 4.4 

3.6 X 108 4.4 
26 Fiber B 1.23 0.63 1.5X10-* 5.0 9.8 

5.1 10° 2.2 
27 Fiber B 2.80 0.63 3.92 x 10-8 5.0 222 

4.3 108 a 
Average Fiber B 

* All experiments are at room temperature. 

rate would be smaller than in the former, due to After loading the bombs, the air was pumped 


a smaller catalytic area per unit volume of out and normal CO, was let in at the proper 
solution. pressure. Samples of the gas were drawn out at 
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ISOTOPIC EXCHANGE 


appropriate intervals (generally from 3 to 2 
half-conversion periods) and analyzed on a 
Nier type mass spectrometer. The constant 
kcC/A was then calculated from (10). Table I 
shows some sample data used in the calculations. 
The equilibrium value of 


(—) fe 
CO,'* normal CO,"* 


was found by leaving the CO, in contact with the 
solution for a long time. Table II gives a sum- 
mary of the results. The precision was limited 
mainly by the spectrometer analyses, but was 
generally within 1'5 percent. 

The interpretation of kcC/A in terms of the 
rate constants in reactions (1) and (2) is made 
by following each isotope separately. For carbon 
exchange : 











k 
C20,4+0H- 2 HC"0,;-, (1) 
Ry 
ky! 
CO: + OH" = HCOs, (1’) 
= 
ke 
CYO.+H:0 2 H.C"0;, (2) 
~2 
ke! 
C*Or+H.0 = H.C"0s, (2') 
hikes ko! k_» 
a= >, ae= . (1 1) 
hb, koko! 
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Considering only (1) and (1’): 


d[HC®O;-] 
— —=k,[C"O2][OH-] 
t 


—k_,{[HC"O;-], (12) 





d[HC"0;-] 
" 1 , — =k,'[ CO, |[OH- ] 
t 


—k_;’THC*O;-]. (12’) 





If W is the volume ratio of the reaction vessel 
to the solution, 


Hi, [HCO;- ] 











(13) 
A W 
Assuming 
[HC®0;-] [H2C*0O3] 
os = , 
[HCO;-] [He2COz;] 
and with 
[C8O, ] 
1 = ——_—__ 
[COs ] 
dN dl HC*O;-]—NdlLHCO;- ] 14) 
caine ocean » 4 
dt [HCO;- |di 


From (11), (12), (13), and (14): 








HdN 1 rOH-JECO i a wy i. 
— =>— a nin r= 
Ad W . Sigg 
(1—N) ky 
= N(1=m) bt - -*) 
(1 —n) k_,’ 


i * 


which is in nearly the same form as Cohen’s 


TABLE IV. Rate constants for exchange reactions. 











Observer Temp. °C ki -mole -min. k_1-min. k2[H20] -min. k-2-min. 
Faurholt (ref. 3) 0 (5.58 X 104) (3.54 10~*) (0.180) (137) 
18 (24.0 10*) (41.6X 10-4) (1.52) (980) 
25 [42.0 104] [110.0 10-4] [3.50] [2120] 
Brinkman et al. (ref. 4) 0 (5.05 X 10°) (0.156) (102) 
18 (12.3 X 10*) (1.44) (666) 
25 [17.5 x 104] [3.46 ] [1400 ] 
Mills and Urey (ref. 2) 0 (0.162) 
0 0.123 
vt 25 1.65 
his paper : 25 
with Fiber A (2.0X 108) (4.4) (6.2) (5500) 
with Fiber B (4.3 X 108) (9.5) (3.3) (2900) 











All values in parentheses ( ) are measured in buffered solutions. 
Values in brackets [ ] are extrapolated. 



































REID AND H.C. UREY 
TABLE V.* 
Diameter Height Flow r(z) ¥(3) 
Pres- Oo ot moles/ r(0) v(0) kcC/A 
Exp. sure Base column column cm? anal- extr. moles/cm* 
no. Packing atmos. used (in.) (ft.) min. ysis tol=% kcCe/L min. 
1 Morell Char. no.7 = 5 2.3 M KOH 0.6 10 1.02 10-3 1.10 1.10 8.0 2.7X10- 
2 Morell Char.no.7 54 2.3 M KOH 0.6 10 1.29x 10-3 1.038 1.038 $4 1.34 1075 
3. 2” Pyrex helices 1 1.3 1 NH,OH 0.6 12 1.38 1.003 1.003 0.096 x 10~° 
4 i”’ Pyrex helices 1.3 2.3 MKOH 0.6 10 0.78 1.11 1.11 8.8 2.3X 10-5 
5 etched 5 4 2.3 M4 KOH 0.6 10 0.78 1.14 1.14 11 2.8X 1075 
6 min,** 50 20% ethan- 
domine 0.6 10 0.98 1.064 1.064 ae 1.7107 
7 i”’ Pyrex helices 1.8.1 NH,OH 0.6 11.8 1.06 1.093 1.103 8.2 2.4K 10-5 
8 etched 23 min. 2.0 1. NH,OH 0.6 12 4.00 1.028 1.035 2.9 3.2107 
9 2” Pyrex helices 
etched 23 min. 1 5.1 14 NH,OH 0.6 11.8 4.00 1.032 1.035 2.9 3.2107 
10 3” Pyrex helices 
etched 23 min. 18 2.6 M NH,OH 0.6 12 8.50 1.020 1.027 2.2 5.1107 
11 3” Pyrex helices 
etched 23 min. 18 2.6 11 NH,OH 0.6 12 17.3 1.044 1.100 8.0 3810-5 
12 3” Pyrex helices 
etched 23 min. 50 2.5 11 NH,OH 0.6 12 16.4 1.083 1.1100 8.8 40x10"™ 
13. 3” Pyrex helices 
etched 23 min. 50 2.5 MNH,OH 0.4 3 34.2 1.012 .013 1.1 40x 10-5 
14. 6-8mesh Haydite 16 2.4 4M NH,OH 0.6 12 15.9 1.032 1.032 2.64 11.5107 
15 Haydite6-8mesh 16 2.4 M NH,OHt 0.6 12 15.9 1.032 1.032 2.64 11.510 
16 Etched 6-8 mesh 
Haydite 16 2.3 M NH,OH 0.6 12 16.1 1.074 1.100 8.0 35.4X 10 
17. Etched 6-8 mesh 
Haydite 6 2.3 MNH,OH _ 0.6 12 15.4 1.030 1.030 23 10.5 1075 
18 Etched 6-8 mesh 
Haydite 33 2.4 WM NH,OH 0.6 12 iss 1.093 1.100 8.0 33.8 10 
19 Etched 6-8 mesh 
Haydite 16 2.3 MNH;OH _ 0.6 1.055 1.055 4.5 8.7X 10 











* All runs are at room temperature. , 
** All etching done with 8.4 per cent HF solution. 
t+ Agitated with etchings. 


equation. When k_, and n and WN are small: 


HH, oN ky 
Foy ltt UcOsMarn—N). (16) 
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A 





A similar equation holds for Eqs. (2) and (2’): 


H.dN ke 
~ [H:O][CO.](am—N). (17) 
Aa W 








Assuming that a; and a2 are near enough to 
the same value so that we may use a common 
value a, then 


IT aN 


k,[OH-]+&2[H 9) 
A dt — ( 1 2 2 





(an—N). 
(18) 


A similar equation holds for oxygen exchange. 
We can now equate kcC/A to 


[CO2] 





ky OH- ko Hoo " 
(kiLOH~]+ho[H20 J) ~ 





Then, with our values of kcC/A we can determine 
the rate constants k; and k2[ HO] if we have at 
least two different sets of conditions for a single 
packing. The calculations of those constants are 
shown for the two Pyrex fibers in Table III, and 
a comparison with other values is shown in Table 
IV. The catalysis of reaction (1) is quite large. 

In the above calculations the pH of the water 
solution was taken from Moore’s tables.!° The 
solubility of CO. was taken from Sander’s 
values.'! The constants used were: 


—log Ka 25° = 6.345 —0.119w, (19)” 
where w is the ionic strength. Thus, with w=0, 
[H+ ]LHCO;-] 
[H2CO; }*LCOsz (dissolved) ] 
10M. B. Moore, with J. H. Buchanan, lowa State Col- 
lege J. Sci. 4, 431-40 (1930). 
11 W. Sander, Zeits. f. physik. Chemie 78, 513-49 (1912). 


12T.), A. MacInnes, Principles of Electrochemistry (Rein- 
hold, 1939). 


=4.52X10-' 





a 25° 

















ISOTOPIC EXCHANGE 


and . 
[ COz (dissolved) | 


onainiansine 880. (20)8-18 
[H2COs] 


Column Experiments 


The results from the reaction vessel experi- 
ments showed enough promise so that actual 
fractionation in packed columns was tried. The 
choice of the packing was indicated by the static 
experiments. Activated charcoal was tried; then 
Pyrex helices. Next, the helices were etched so as 
more closely to simulate the Pyrex fiber surface. 
Then a cheap, siliceous, shale packing, Haydite, 
was tried, both unetched and etched. A control 
using the unetched Haydite and a solution 
agitated with etchings was run to make sure that 
the increased rates were not due to homogeneous 
catalysis from some material left on the etched 
surface. The runs were made at room tempera- 
ture. Pressures from one to fifty atmospheres 
were tried. The procedure was straightforward: 
Columns were filled with the packing to be inves- 
tigated. The column was filled with CO, and a 
solution of base was introduced at the top at the 
rate desired, reacting with the CO. to form 
HCO;-, which traveled down the column. Nor- 
mal CO, was fed in at the bottom at constant 


13 E, R. Strohecker, Zeits. f. Nahr.-Genussm. 31, 121-60 
(1926). 

4Q. Haechnel, Centr. Mineral Geol. 1920, 25-32 (1920). 

6 F. J. J. Buytendyke, R. Brinkman, and H. W. Mook, 
Biochem. J. 21, 576-84 (1927). 

1% E. Miiller and A. Luber, Zeits. f. anorg. allgem. 
Chemie 190, 427 (1930). 
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pressure. In these experiments, since the heavier 
carbon concentrates in the liquid phase, the C” 
concentrated at the top. A steady drain took 
place at the bottom while samples of the gas 
were drawn off at the top and analyzed. kcC/A 
was then calculated from (7). A summary is 
given in Table V, illustrated in part by Figs. 1, 
2, and 3. 

In most instances the equilibrium value of 
r(z)/r(0), where r is the ratio of the heavy to the 
light isotope, was determined by analyzing one 
sample near the first part of the run and getting 
checks on the same value from two successive 
samples after equilibrium should have been 
reached. 

Experiment No. 4, a little more detailed, is 
outlined in Fig. 4. The measured holdup of 
solution was 50 cc. From the data given and (9), 
7=4.4 hours; whereas, from the experimental 
curve, 7=5.2 hours, in fair agreement with 
Cohen’s theory. 

Fluid flow theory states that in flow through 
a packed column, a fluid maintains laminar flow 
only up to a given critical value of the Reynolds 
number. When uy is the viscosity coefficient, p is 
the density, U is the linear velocity, and d is a 
measure of the spacing between the packing, the 
Reynolds number is given by 


R.=dU p/p. 


The critical value is around 40, depending on the 
type of packing; however, it is rather loosely 
determined for any particular application. In the 
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Fic. 2. Pyrex helices, etched 23 min., 1 atmos. —o—. 
6-8 mesh Haydite, etched 23 min., 16 atmos. --- X---. 


four cases most widely investigated, shown in 
Figs. 3 and 4, the approximate gas flow rate to 
Reynolds number correspondence was deter- 
mined. Using a value of 0.7 cm ford, R-=3 XK 108 U 
min. cm?/mole. This would place the large in- 
crease of kcC/A per unit flow at a Reynolds 
number from 25 to 50 for all the cases except 
the atmospheric run. From this it seems likely 
that the establishment of turbulence is respon- 
sible for the increase, and since much slower 
rates of flow have produced faster exchange in 
other systems, it appears that diffusion through 
the liquid to the catalytic surface is rate con- 
trolling, and may be speeded up by producing 
turbulence in the gas. 

From the results obtained, the specifications 
of an efficient system for fractionating carbon or 
oxygen isotopes by the above method should 
include: 
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Fic. 3. Pyrex helices, etched 23 min., 18 atmos. —o—. 
Pyrex helices, etched 23 min., 50 atmos. --- X---. 





I. A pressure somewhat above 6 atmospheres 
up to 18 atmospheres so that [COz (dis.) ] 
is great enough. 

II. A high enough fH so that advantage may 

be taken of catalyzed (1). 

An economical packing (such as etched 

Haydite) with efficient fractionating charac- 

teristics and sufficient catalytic activity so 

that (1) and (2) are not rate controlling. 

A gas flow of sufficient velocity for turbu- 

lence. 

V. As high a temperature as is necessary to 
keep the rate of (3) at the required value, 
but not too high, as @ decreases with 
increasing temperature.’ 


Il. 


IV. 


Extension of Theory and Results to 
Large Scale Operation 


Using the above results we calculated the 
specifications for erection and operation of a 
column system which would produce about 4 
grams of carbon, 33 percent C™, per day. In 
order to concentrate C" to where it is about 50 
percent of the C™ concentration, with a=1.012 
and assuming normal COs, to be 1.10 percent C™, 
320 plates would be necessary. Because of the 
time required for a large single-stage system to 
come to equilibrium, usually a cascade system 
is designed for such high concentration. 

This calculation is on such a six-stage cascade 
with the fractionation split evenly between the 
stages. The design is such that each unit has the 














es 


)] 


mM 


de 


he 











° 
° 
o- 
— a 
° 
co} y 
- i 
> oO 
N 
. 
o 
oL 
° 
° 1 1 — 
0 10 





20 Operation Time far) 


Fic. 4. Column experiment No. 4. 0.78 mole/cm? min. flow. 


same C" transport with a material flow such 
that 80 percent of the maximum C* possible is 
stripped. The height of each column is calculated 
by (8). 

The calculation of the time required for the 
cascade to come to production cannot be made 
by Cohen’s theory, so another fairly accurate 
approximation was made. This consisted of de- 
termining the amount C™ over the normal there 
was in the cascade and dividing it by the rate 
at which C® was stripped in the first column. If 
Q is the total moles of C® in a column approxi- 
mately 


O=(h+H) f Na: (21) 


N r.\%l8 
r=( ) and v= (“) , 
1—N ro 


assuming r,=frov?, (21) becomes 


With 





z Yo 72 
Q=(h+H) “ype 


0. ro 


sian (—") 
= n e 
In V 1+79 





ve 


(22) 





If R is the number of moles of C™ stripped off in 
the first column A, then 


T 
R = Lif (N. nas n2a)dt, (23) 
0 


where T is the time required for equilibrium. If 


i N.a(1 —MN2A) 


~ ma(1— Nea) 
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TABLE VI. Specifications for a 6-stage cascade. 16 atmos- 
pheres COz with etched Haydite packing. 





























Inside Flow of % Cw 
diam. Height CO2 at 
Unit (inches) (feet) moles/min. bottom 
A 3.00 78 0.726 2.05 
B vB «| 82 0.395 3.80 
c. 1.64 85 0.218 6.93 
D 1.23 87 0.123 12.3 
E 0.95 89 0.073 20.9 
F 0.77 91 0.048 33.0 
TABLE VII. 
Cyanide Cascade Bicarbonate Cascade 
1.012 1.012 
8.11074 3.51074 
kcC/A moles cm=3 moles cm 
(operating) min.—!* min.~} 


78’X 3" diam. 

82’ 2.21" diam. 
85’ 1.64” diam. 
87'X 1.23” diam. 
Size of 5th unit 89’ 0.95”" diam. 
Size of 6th unit 91’X0.77” diam. 
Flow through 1st 0.042 mole min. 0.726 mole min.~! 

unit min.~! 
Production of C® 0.161 g per day 
13 


40’ 1" diam. 
25’ xX 2” diam. 


Size of 1st unit 
Size of 2nd unit 
Size of 3rd unit 
Size of 4th unit 


1.46 g/day 

ete 0.23 0.33 

Time to reach 12 days 87 days 
production 





* Value for the 1” column. For the 3” column, kcC/A =2.8 X1073 
mole cm=* min.~. 





and N.,4 is small 


TB-1 
R ae Dal os ——dt, (24) 
0 8B 
and assuming a mean value for B: 
(8-1) 
R=LsaN.a——T. (25) 
B 


We have for the complete system: 
2Q0;:=R+Near(hitH,)Zi, (26) 


which becomes 
B 
T =————22,(h; + Hi) 


(8@—1)La 
i 1— No; 
In (—*) 
1-—N,; 


x =—l. @n 





N2i(1— Noi) 
N.a ln ——————_ 
a Noi(1—N3:) J 
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In the above system using 6-8 mesh Haydite 
in the larger columns and 16 atmosphere of 


COs (cf. Table VI) 


IT/A =0.84X 107 mole/cm’, 
h/A =0.42X 10 mole/cm*. 


The flow planned was 0.0160 mole/cm? min. 
Under these conditions a value of 3.5X10- 
mole/cm* min. was used for kcC/A. 

With p4=0.000328, the daily production is 
1.46 g of C. The time to come to production, 
from (27), assuming B= 1.011, is 87 days. 


HUGGINS 


Comparison of Cyanide and 
Bicarbonate Methods 


A comparison of the HCN—CN-'® and 
CO2z—HCO;- methods of C® concentration can 
be made by a glance at Table VII. For large scale 
operation, the cyanide method has the advantage 
of a smaller holdup and thus a shorter length of 
time to reach production. On the other hand, the 
raw materials are more expensive, and it has the 
aforementioned undesirable properties of poly- 
merization and toxicity. In either case, the time 
necessary to reach production could be reduced 
by using a larger number of units in the cascade. 
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The Photographic Latent Image 


I. Lattice Energy and Related Quantities for a Hypothetical Silver Bromide 
Crystal Having the Cubic ZnS Type of Structure*; 


Maurice L. HuGGins 
Eastman Kodak Company, Rochester, New York 
(Received April 2, 1941) 


The usual NaCl type of structure (B/) for silver bromide is compared with a hypothetical 
structure for this substance of the cubic ZnS type (B3) with regard to certain quantities which 
are of importance in connection with the theory of the photographic latent image. These include: 
lattice energy, equilibrium distance between closest atoms, density, potential at interstitial! 
positions, potential energy changes as an ion moves from its normal position to an interstitial 
position, and as an ion or electron moves from one interstitial position to another. The methods 
of calculation are outlined and the probable accuracy of the results discussed. 


INTRODUCTION 


HE second paper of this series will contain a 

critical discussion of the hypothesis that, 
in the photographic latent image in a silver 
bromide crystal, the region immediately sur- 
rounding the photoelectrons has a structure. of 
the B3 (cubic ZnS) type—or another structure 
closely related to it—instead of the BJ (NaCl) 
type structure possessed by the rest of the 
crystal (see Fig. 1). 


* Presented before the Division of Physical and Inorganic 
Chemistry at the St. Louis Meeting of the American 
Chemical Society, April 9, 1941. 

+ Communication No. 802 from the Kodak Research 
Laboratories. 





Preliminary to this discussion it seems ad- 
visable to outline the methods and present the 
results of certain pertinent calculations, e.g., cal- 
culations of the approximate energy and dimen- 
sional changes to be expected from the B1 to 
B3 transformation and of the energy changes to 
be expected during the passage of Ag+ ions, Br- 
ions, or electrons through each of these types of 
structure. 


LATTICE ENERGY OF AgBr HAVING 
THE Bl TYPE STRUCTURE 


Mayer! has computed a theoretical value of 
13.5-10-' erg per molecule or 197 kcal. per mole 


1J. E. Mayer, J. Chem. Phys. 1, 327 (1933). 
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for the lattice energy, to be compared with an 
‘experimental’ value, computed by means of a 
“Born cycle,’’ of 202 kcal. per mole. He suggests 
that the difference is primarily the result of the 
neglect in the theoretical calculation of any 
homopolar bond energy contribution. 

In Mayer’s treatment, certain constants in the 
exponential expression for the “‘overlap’’ repul- 
sion energy between two ions were derived from 
Richards and Jones” compressibility data for 
silver halide crystals. In view of Bridgman’s* 
more recent results on these crystals, it seems 
that the probable accuracy of the data hardly 
warrants this procedure, rather than the simpler 
alternative of assuming the constant a, in the 
expression , 


exp La(ragtrsr—?) ], 


to be the same as that (3.00-108 cm~') found to 
be satisfactory for all the alkali halides. This 
latter procedure, moreover, has the advantage 
of giving a value of ra, which is directly com- 
parable with the other basic radii previously 
obtained.* ® 

A recalculation of the lattice energy of silver 
bromide (B1 type) has, therefore, been made, 
following strictly the procedure used® for the 


*T. W. Richards and G. Jones, J. Am. Chem. Soc. 31, 
158 (1909). 

*P. W. Bridgman, Proc. Am. Acad. 74, 21 (1940). 
(193) L. Huggins and J. E. Mayer, J. Chem. Phys. 1, 643 
5M. L. Huggins, J. Chem. Phys.’ 5, 143 (1937). 
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Fic. 1. Unit cubes of the B/ and B3 structures. 


alkali halides and adopting Mayer’s estimates of 
the van der Waals attraction constants. The 
value obtained is 13.2-10-” erg per ‘‘molecule.”’ 
The difference between this and Mayer’s value 
is certainly within the probable error of the cal- 
culation. 

The value of ra, deduced in this calculation is 
1.05-10-§ cm. Using the approximate equation 
previously given [Eq. (5) of reference 5], with 
appropriate values of the constants, the com- 
pressibility is computed to be 2.4-10-" barye™. 
This is in sufficiently good agreement with 
Richards and Jones’ value of 2.74-10-" and 
Bridgman’s value (which, he states, is probably 
too low) of 2.30-10-". 


NON-COULOMB INTERACTION ENERGIES 


If the same constants are used as in the lattice 
energy calculations, the non-Coulomb interac- 
tion energy between two Ag* ions varies with 
the distance r between them according to the 
equation 


bagdg= — 67 -10-%7-*— 91 - 10-78 

+1.111 exp [3(2.10—10*7) ]. (1) 
Similarly, 
dagnr= — 109+ 10-8%-8— 199. 10-78y-8 

+0.9305 exp [3(2.61—108r)] (2) 
and 


dprpr= — 208-10-*r-*— 475 - 10-8 
+0.750 exp [3(3.12—10*r)]. (3) 








414 


Curves representing these three equations are 
shown in Fig. 2. Within the error’of the calcu- 
lation, @agag is zero for all interionic distances 
greater than 2.5A. 

It should be noted that these equations and 
curves make no allowance for the polarization of 
each ion by the charge on the other. This polar- 
ization would increase the repulsion between 
two like-charged ions and increase the attraction 
(decrease ¢) between two oppositely charged 
ions. For structures in which each positive 
ion is symmetrically surrounded by negative 
ions and vice versa, this polarization can prob- 
ably be neglected. 


LATTICE ENERGY AND INTERATOMIC DISTANCES 
FOR AgBr HAVING THE B3 TYPE STRUCTURE 


Silver iodide at room temperature exists in 
two closely related types of structure, desig- 


107"“erg 




















Fic. 2. Non-Coulomb interaction energies between pairs of 
ions, neglecting polarization. 
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HUGGINS 


nated® as B3 and B4, of nearly equal stability ; 
in both types, each ion is tetrahedrally sur- 
rounded by four of the other kind. It is of interest, 
for our present purpose, to compute the lattice 
energy of a hypothetical silver bromide crystal 
having the B3 structure, for comparison with 
the lattice energy given above for the Bi type. 

The lattice energy depends on the scale of the 
structure, hence on the equilibrium distance be- 
tween closest atoms. This distance can be com- 
puted, approximately, by the method of Huggins 
and Mayer [reference 4, Eq. (13) ], using the 
same constants for the interactions between 
pairs of ions as in the calculations for the Bl 
type. The assumption is involved that the struc- 
ture can be treated as if purely ionic. This is 
doubtless a poorer assumption than for the B/ 
type and some error is introduced on that ac- 
count. The true equilibrium distance should be 
less than that computed. 

The ratio of the thermal expansion coefficient 
to the compressibility also enters into the calcu- 
lation, but large differences in the assumed ratio 
do not greatly affect the results. If one takes this 
ratio to be the same as for the BJ type (expan- 
sion coefficient from International Critical Tables, 
compressibility from Richards and Jones’), the 
computed value of 7 (the equilibrium distance 
between closest ionic centers) is 2.73A (and so 
ao= 6.305A). 

The corresponding molecular volume and 
density are 62.7A* and 4.98 g/cc, respectively. 
These are 13 percent greater and less, respec- 
tively, than the corresponding quantities for the 
B1 type. The computed value of the lattice en- 
ergy for the B3 structure is 12.79-10-. Com- 
parison with the calculated lattice energy for the 
B1 structure shows a difference in stability of 
about 0.4-10-” erg per molecule, i.e., about 0.3 
electron-volt per molecule or 6 kcal. per mole. 
Since the B3 structure (but not B/) is such that 
each atom can be bonded at the same time by 
(partially polar) covalent bonds to each of its 
close neighbors, one would expect the neglected 
contribution of the covalent bond energy to 
make the lattice energies—and so the stabilities 
—of the two types even more nearly equal. 


6P, P. Ewald and C. Hermann, Strukturbericht, 1913- 
1928 (Akademische Verlagsgesellschaft, Leipzig, 1931). 








THE 


POTENTIAL ENERGY CHANGES FOR MOTION OF 
IONS AND ELECTRONS IN THE BI STRUCTURE 


It is desired to compare the potential energy 
changes in the Bl and B3 structures, as elec- 
trons or Ag+ or Br~ ions move along the lines of 
least resistance through them. These energy 
changes can be approximated by considering the 
electronic or ionic motions to occur without dis- 
placement or polarization of the other ions. This 
undoubtedly introduces considerable error, but 
the errors should be in the same direction and 
roughly of the same magnitude (or the same frac- 
tion of the total energy), at corresponding points 
in the two structures. 

The ionic motions considered are those in- 
volved in the production of interstitial ions and 
in the motion of such ions from one interstitial 
position to another. Estimates of the net energy 
change when an ion moves from a normal lattice 
position to a distant interstitial position and of 
the magnitude of the energy hump between two 
interstitial positions have previously been made* ® 
for NaCl, KCI, and RbCl. Numerical comparison 
of these estimates with those deduced from the 
work reported here for AgBr will not be pre- 
sented, because of the differences in the sub- 
stances, because of important differences in the 
expressions used for the computation of the non- 
Coulomb energies, and because of the approxi- 
mations, in the AgBr calculations, mentioned in 
the preceding paragraph. 

The energies of interest are in part due to 
Coulombic and in part to non-Coulombic inter- 
actions. For the non-Coulomb part, we use 
Eqs. (1) to (3), summing over all the ions closest 
to the ion in question. To obtain the Coulombic 
contribution, we compute the Coulomb poten- 
tial at appropriate points along the lines of 
motion in which we are interested and then 
multiply by the charge (+e) on the moving 
particle. 

The treatment of an electron as a point charge, 
moving along certain fixed paths in the structure, 
is doubtless quantitatively inaccurate, but should 
give approximately correct relative values for 


on Jost and G. Nehlep, Zeits. f. physik. Chemie B32, 1 


(1 
*N. F. Mott and M. J. Littleton, Trans. Faraday Soc. 
34, 485 (1938). 
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IMAGE I 


Fic. 3. Illustrating H¢jendahl’s method of computing the 
potential at a point in a B/ structure. 


the energy changes accompanying motion in the 
two crystal types being discussed or along differ- 
ent paths in either. 

For the BI structure, we are interested in the 
motion of an Ag+ or Br~ ion from its normal 
position along a cube diagonal to the center of an 
adjacent eighth-unit cube, also in the motion of 
an ion or electron from the center of one such 
small cube through one of its faces to the center 
of an adjacent cube (see Fig. 1). 

The Coulomb potential, due to the other ions, 
at the center of any one ion in its normal posi- 
tion is tae/ay, where a is Madelung’s constant 
(=3.4951 for Bl), e is the magnitude of the 
electron charge, and a» is the unit distance 
(=2ro for B1). The Coulomb potential at any 
point on the straight line connecting the centers 
of two adjacent cubes is zero, if no ions are miss- 
ing or displaced from their normal positions. 
In case an ion is missing from the structure, the 
potential at a point on this line is +e/d, d being 
the distance between the vacant position and the 
point in question. 

To obtain the Coulomb potential at a point 
on the diagonal connecting a normally occupied 
position with a cube center, we follow the method 
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Fic. 4. Calculated energies for motion of an ion or 
electron in a B/ type AgBr crystal along certain paths con- 
necting points a, d, b, and c in the lattice (see Fig. 1). The 
ordinate scale at the upper right applies only to curves VJ 
and III (between 6 and c). For further explanation, see 
Table III and the text. 


used by Hgjendahl® in computing the Madelung 
constant. We sum 1/d directly, collecting to- 
gether all the terms for ions at the corners of the 
eighth-unit cube containing the point for which 
the calculation is being made, then collecting 
all the terms for ions in the faces, edges, or 
corners of the next cubic box, concentric with the 
first, etc. (see Fig. 3). Three such boxes suffice 
to give an accuracy of better than 0.03. 

Multiplying the Coulomb potentials, obtained 
as just described, by the magnitude of the elec- 
tron charge gives the potential energy. This 
changes along the lines of interest as shown by 
curve I of Fig. 4. 


°K. Hgjendahl, Kgl. Danske Vid. Sels. Math.-Fys. 
Medd. 16, 133 (1938). 




































HUGGINS 


For the motion of a silver or bromide ion from 
a corner to the center of an eighth-unit cube or 
from one cube center to an adjacent one, the 
total energy—obtained by adding, to the Cou- 
lomb energy, terms for the non-Coulomb inter- 
action with the closest other ions—varies as 
shown by the other lines in this figure. 


POTENTIAL ENERGY CHANGES FOR MOTION OF 
IONS AND ELECTRONS IN THE B3 STRUCTURE 


The motions of interest in the B3 structure are 
all along the body diagonals of the cubes. Be- 
sides the positions (a’ and d’, Fig. 1) normally 
occupied by ions, two types of normally unoc- 
cupied positions are of especial interest: (1) the 
b’ positions, which are tetrahedrally surtounded 
by four Ag+ ions, each at the distance ro (the 
same as the distance between closest ion centers), 
and (2) the c’ positions, tetrahedrally surrounded 


by four Br- ions, each at the distance 7. The’ 
J ’ 


b’ positions are at the centers of eighth-unit 
cubes, if the unit cell is so chosen that an Ag* 
ion is at each cell corner, or at the center and 
the centers of the edges of the unit cell, if it 
is so chosen that a Br~ ion is at each cell corner. 
Likewise, the c’ positions are either at the center 
and the centers of the edges of the unit cell or 
else at the centers of half of the eighth-unit 
cubes, depending on the arbitrary choice of 
origin of the unit. 

We are interested in the motion of an ion 
from its normal position (e.g., a’) along the cube 
diagonal to the center (b’) of an adjacent eighth- 
unit cube and then to the center (c’) of the unit 
cell. We are also interested in the motion of ions 
and electrons along the cube diagonals (b’ to c’) 
in a B3 structure in which no ions are missing 
from their normal positions. 

The Coulomb potentials at points along the 
diagonals cannot readily be calculated by Hgjen- 
dahl’s method. They have been computed, how- 
ever, by a modification of Madelung’s method."” 
We first choose an orthorhombic unit cell for the 
structure, with a diagonal of the cubic unit cell 
as the X axis of the new unit (see Fig. 5). The 
coérdinates of the ions are listed in Table I. 

We wish to calculate the potential at a point, 
% 00, with x) having values between 0 and 


10 E. Madelung, Physik. Zeits. 19, 524 (1918). 
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TABLE I. Coordinates of ionic centers, referred to the 
orthorhombic unit cell of Fig. 6. 





| Ag* Br- 
x y z x y 4 
0 0 0 1/4 0 0 
2/3 1/2 1/6 11/12 1/2 1/6 
1/3 0 1/3 7/12 0 1/3 
¢ 0 1/2 1/2 1/4 1/2 1/2 
lols 2/3 0 2/3 11/12 0 2/3 
" 1/3 1/2 5/6 7/12 1/2 5/6 


TABLE II. Values of x, y and z for substitution in 











y x Eqs. (6), (7), and (8). 
4%.0 -— a=5a, oy - so ee 
i ie — N ; P 
Fic. 5. Orthorhombic unit cell for the B3 structure, used Ee EL 
in calculation of the changes in Coulomb potentialasanion 4, 7, 13 2/344 1/2 1/6 
moves along a cube diagonal (=the X axis). 2, 8, 14 1/3+x¢ 0 1/3 
3,9, 15 Xo 1/2 1/2 
i ‘ pee One 4, 10, 16 2/3+X0 0) 2/3 
1/2. (The ion at the origin is assumed to have | 5) 44) 17 1/3+%o 1/2 5/6 
been removed.) 6, 12, 18 Xo 0 I 
The Coulomb potential due to other ions in) === ees —omeneenene 
the X axis is given by the equation 
e 1 10~* ere 10°'*erq 
Vo ie ——— ; 1s 7's 
a xot} 


vi 





x 1 1 
m=1 16m—16x»9—4 16m—16xo 


1 1 
creme reser seenatere -)} (4) 
16m+16x, 16m+16x9+4 
Calculations for values of x» which are multiples 
of 1/16 are simple, since the denominators of the 
terms to be summed are all integers, and tables 
of reciprocals of integers are readily available. 
We next consider the Coulomb potential at 
t00 due to other rows of. ions, parallel to the 
X axis, in the X Y plane. This is given, for the 
nth row on either side of the X axis, by the 
equation 


4e , 
V,=— > Ko(2aIrb/a) 
a l=1 


- feos (2alxo) —cos [2l(xo+})]}. (5) 














The function Ko (x) is that tabulated by 
. ° . ‘tor d’ ’ : ' 
Jahnke and Emde" as (ix/2)Ho (ix). A very ow ee os 
small number of terms suffices. Fic. 6. Calculated energies for motion of an ion or 
Sa electron in a B3 type AgBr crystal along certain paths 
"E. Jahnke and F. Emde, Funktionentafeln mit Formeln connecting points a’, d’, b’, and c’ in the lattice (see Fig. 1). 
und Kurven (Teubner, Leipzig and Berlin, 1909), p. 135. For further explanation, see Table III and the text. 
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TABLE III. Summary of energy changes (in 10-" erg), computed for undistorted structures and neglecting polarization 











Curves, Net energy Activation 
Figs. 4 increase energy 
and 6 1 B3 B3 
Motion of Agt from normal position to nearest stable interstitial rej 10.9 8.3 13.2 8.5 
position. 
Motion of Br~ from normal position to nearest stable interstitial III 26 17 42 2Z 
position. 
Motion of electron from one (most stable) interstitial position to IV 0 0 0 2.1 
another like it, in a complete lattice. . 
Motion of Ag* from one (most stable) interstitial position to another V 0 0 3.6 is 
like it, in complete lattice. 
Motion of Br~ from one (most stable) interstitial position to another VI 0 0 18 7 
like it, in a complete lattice. 
Motion of Agt from one of the interstitial positions of greater stability V -— 3 -- 1.5 
to one of those of lesser stability, in a complete crystal. 
Motion of extra electron from an interstitial position in B/ to one —1.1 
in B3. 
Motion of extra Agt from an interstitial position in BJ to one in B3. —2.8 
Transformation, BI to B3 (per ‘‘molecule’’) <0.4 








The Coulomb potential due to ions in the 
Nth plane, parallel to the XY plane, on each 
side of that through the origin, is 


2-6'e 
Vy =——{2 © & Kim cos 2xmy 


a l m 


X [cos 2xlx —cos 2nl(x+4) ] 
+> Kif[cos 2rlx—cos 2l(x+4)]}, (6) 
l : 








where 
exp [ —2)22(/?+6m?)? ] 7) 
‘in (2-+6m?)! 
exp [ — 2x21] 
K, = ] ’ (8) 


and x, y, and z have the values listed in Table IT. 
The total Coulomb potential due to all the 
ions in the lattice (other than that at 000) is 


V’=Votd Vatd Vn. 


We have computed it for values of xo from 0 to 
1/2, in multiples of 1/16. For x»=0, we obtain 
V’ = —3.7829e/ao, the constant agreeing with the 
value of Madelung’s constant for this structure.” 

The addition of e/xoa to the calculated values 
of V’ gives the Coulomb potentials V at the 
points in question for the complete structure 
(including a positive ion at 000). As should be 
the case, V is symmetrical about x»= 3/8, where 
it is zero. (See curve IV, Fig. 6.) 





12Q, Emersleben, Physik. Zeits. 24, 97 (1923). 


At x9=0.25, V has a positive value of 0.500e/a 
(=0.2887e/ao). A crystal having a9=6.305A, 
with an electron in such a position, is 1.06-10-” 
erg more stable than if the electron were not 
there or were in a position of zero potential. 
At x9=0.50, V has the same magnitude, but is 
negative. 

We now add to the Coulomb energy the sum 
of the non-Coulomb interaction energies for all 
the closest atoms (Fig. 2) to get the total energy 
for a silver or bromide ion at any given position 
on the cube diagonal. In Fig. 6 are curves show- 
ing how the total energy changes as an ion 
moves from its normal position (a’) along a 
diagonal to nearby interstitial positions (b’ and 
c’), also similar curves for motion of an extra 
ion or electron from one interstitial position 
(b’ or c’) to another, in a complete lattice. 
These curves are all computed for the interionic 
distance, 7o=2.73A (or a9=6.305A). 


SUMMARY AND CONCLUSIONS 


The results of the calculations reported here 
are summarized in Figs. 4 and 6 and in Table III. 

Corrections for polarization and distortion of 
the structure in the vicinity of extra ions, extra 
electrons, and vacant normally-occupied posi- 
tions would probably decrease all the net energy 
increases and activation energies in this table. 
The approximate relationship between the Bl 
and B3 values should not be markedly affected, 
however. 

It should also be remarked that the energy 






























THE 





changes indicated for the motion of an electron 
through either structure do not include any 
“overlap” energy of repulsion between the elec- 
tron and neighboring atoms, nor is any considera- 
tion given to possible discrete energy levels for 
this electron. An electron moving through a Bl 
type structure ‘passes much closer to other atoms 
than one moving through a B3 structure; this 
might result in a considerably larger ‘‘activation 
energy’ hump in the former case than in the 
latter. 

From the results of these calculations, it seems 
safe to draw the following conclusions: 
1. For a perfect crystal, the B3 structure is 


PHOTOGRAPHIC 





LATENT IMAGE II 419 








slightly less stable than the B/ structure. This 
relationship may be reversed, however, in the 
vicinity of an extra Ag+ ion or an extra electron, 
or both. 

2. Ionic conduction through both structures is 
by Ag* ions, not by Br- ions. 

3. The ionic conductivity of the B3 structure 
is considerably greater than that of the B/ struc- 
ture, at the same temperature.'® 

4. Extra electrons and extra ions ‘“‘prefer”’ 
the B3 structure to the B/ structure. 


13 Tt may be noted that AgI, with a B3 or B4 structure, is 
a much better conductor than AgBr, with a B/ structure, 
at the same temperature. 
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A theory of the formation and growth of the photographic latent image and of the develop- 


ment process in the crystal containing it is outlined. It is assumed that the latent image consists 
of a small region of the (silver halide) crystal in which the atomic arrangement is the ‘‘tetra- 
hedral’’ B3 type (or another closely related arrangement), stabilized by the presence of one or 


INTRODUCTION 


HEN a silver bromide crystal held in a 

gelatin layer is exposed to light, some 
change occurs, as a result of which a developing 
solution can attack the crystal and reduce it to 
metallic silver more easily than if it had not been 
exposed to light. In photographic literature, a 
crystal which has suffered this change is said to 
contain a latent image. Many theories have been 
suggested as to the nature of the change which 
occurs.' That most generally accepted at the 


* Presented before the Division of Physical and Inorganic 
Chemistry at the Detroit Meeting of the American 
Chemical Society, April 12, 1943. 

+t Communication No. 922 from the Kodak Research 
Laboratories. 

*C. E. K. Mees, The Theory of the Photographic Process 
(Macmillan, New York, 1942). 


more photoelectrons, with or without extra silver ions. 





present time is due to Gurney and Mott.?* Ac- 
cording to this theory, the initial effect of light is 
to liberate photoelectrons in the silver bromide 
lattice. These move to specific points and produce 
a negative charge at those points. Silver ions then 
move through the lattice; when they reach the 
charged points they are neutralized, forming a 
metallic silver nucleus. This silver nucleus, it is 
assumed, constitutes the latent image, which 
makes the crystal more readily developable than 
it would be otherwise. 

The purpose of the present paper is to suggest 
an alternative theory, according to which the 
latent image consists not of crystalline metallic 


2R. W. Gurney and N. F. Mott, Proc. Roy. Soc. 164A, 
151 (1938). 
3 J. H. Webb, J. App. Phys. 11, 18 (1940). 
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silver but of a small region of the crystal or 
crystal surface in which the atomic arrangement 
is that known as B3 or cubic ZnS type (or some 
other closely related arrangement), stabilized by 
the presence of extra electrons (photoclectrons), 
perhaps with extra silver ions also. 


LATTICE ENERGIES 


AgCl and AgBr crystals normally have the 
NaCl or B1 type of structure,‘ in which each ion 
is surrounded equidistantly by six ions of the 
other kind. AgI crystallizes in either of the two 
closely related ‘‘tetrahedral’’ structures, B3 
(cubic) and B4 (hexagonal), in each of which each 
ion is surrounded tetrahedrally by four of the 
other kind. (For diagrams of the unit cubes of the 
B1 and B3 structures, see Fig. 1 of the preceding 
paper.®) AgBr and AglI form solid solutions, 
having the B/ structure if the iodide concentration 
is low, the B3 or B4 structure if it is high, and 
either or both structures within an intermediate 
range.® 

The lattice energies of the B3 and B4 types of 
structure are very nearly the same. For simplicity 
in description, we shall consider, in the remainder 
of this paper, only B3 (of these two). It should be 
understood, however, that nearly all that is 
stated with regard to B3 would also apply to B4. 

From the facts just noted about the structures 
of AgBr and AglI and their solid solutions, one 
would expect the lattice energies of a BJ type 
AgBr crystal and of a hypothetical B3 type AgBr 
crystal to differ but slightly. Calculations of these 
energies,> 7 by methods which have previously 
yielded results in good agreement with experi- 
ment, show this to be the case. The lattice 
energies for the two structures should be even 
more nearly equal for silver bromide containing 
a small amount of iodide, asin many photographic 
emulsions. The preference of the iodide for the 
tetrahedral structures may be explained as due to 
the greater tendency of iodine (than of bromine) 
to form covalent bonds with silver. Electron-pair 
bonds between all adjacent atoms are possible in 
the B3 and B4 structures,* but not in the Bl 


*P. P. Ewald and C. Hermann, Strukturbericht, 1913- 
1928 (Akademische Verlagsgesellschaft, Leipzig, 1931). 
5M. L. Huggins, J. Chem. Phys. 11, 412 (1943). 

*R. B. Wilsey, J. Frank. Inst. 200, 739 (1925). 
7J. E. Mayer, J. Chem. Phys. 1, 327 (1933). 
5M. L. Huggins, J. Am. Chem. Soc. 44, 1841 (1922). 
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structure, the number of valence electrons being 
insufficient. 

The B3 structure, unlike the B/ structure, 
contains “‘holes”’ of positive potential and other 
“holes” of negative potential (see Fig. 1). The 
former are possible sites for extra electrons (e.g., 
photoelectrons). From the calculations of the 
preceding paper, we may conclude that an 
electron in such a site stabilizes the structure by 
at least 1-10-" erg or 0.6 electron volt. The 
presence of extra electrons therefore favors the 
transition from a B1 to a B3 structure. 

The holes of positive (and negative) potential 
in the B3 structure are also possible sites for extra 
negative (and positive) ions. Extra ions (if not 
too large) would fit into the crystal structure 
without introducing large overlap repulsion 
energies. An extra Agt ion, for instance, would 
have as nearest neighbors four tetrahedrally 
disposed Br~ ions, these closest ions being ar- 
ranged just like those around each of the other 
similar normally placed ions in the structure. 

_ Extra ions cannot fit into the BJ structure 
without producing large repulsion energies. From 
the data in the preceding paper, one may com- 
pute 2.9-10-™ and 0.15-10- erg as the increase 
in energy (decrease in stability) resulting from 
the insertion of an extra Agt ion in the B/ and B3 
structures, respectively. Allowance for polariza- 
tion and shifts of the surrounding ions would 
decrease these instability energy values some- 
what; the actual B3 energy change is, in all 
probability, negative, the structure being more 
stable with the extra Agt ion than without it. 
(Many extra Ag* ions cannot be present, without 
compensating extra electrons or negative ions, 
because of the Coulomb energy produced by their 
mutual repulsions.) 

Similar calculations indicate that the chance of 
having extra Br~ ions in either type of structure 
is negligible, except possibly at the surface of a B3 
structure, where large displacements can occur. 
Additional Br~ ions might well be adsorbed at 
otherwise vacant surface sites of positive 
potential. 

One may conclude from the foregoing that 
extra Ag* ions within the crystal or extra Agt 
ions or (possibly, to a lesser extent) Br~ ions on 
the surface tend to make the B3 structure more 











ing 


Ire, 
her 
Phe 
2, 
the 

an 
by 


The 
the 


tial 
‘tra 
not 
ure 
sion 
yuld 
ally 
ar- 
ther 
e. 
ture 
rom 
om- 
ease 
rom 
1 B3 
riza- 
ould 
yme- 
1 all 
nore 
t it. 
hout 
ions, 
their 


ce of 
‘ture 
a B3 
ccur. 
d at 
itive 


that 
Agt 
is on 
more 








= 


PHOTOGRAPHIC 





LATENT IMAGE Il 























4 
4 
4 
‘ 
‘ 
4 















































@ Ag” 
OBr 


e SITE FOR EXTRA + ION 


o SITE FOR EXTRA —10N OR ELECTRON 


Fic. 1. A portion of the B3 structure, showing (small dots) possible sites for 
extra positive ions and (small open circles) possible sites for extra negative ions or 


electrons. 


stable than the B/ structure, in their immediate 
vicinity. 

Extra ions of other sorts—especially those 
having high polarizabilities and tendencies to 
form covalent bonds but with radii about the 
same or smaller than that of Ag+—should simi- 
larly favor the B3 structure. Particularly effec- 
tive should be the substitution of a bivalent 
positive ion, such as Hg**, for some of the Ag* 
ions, along with extra electrons—one for each 
substitution, if the crystal remains electrically 
neutral. The addition of mercury atoms (Hg** 
plus two electrons) to the surface of a silver 
halide crystal would certainly favor the transition 
from BI to B3 in the adjacent positions of the 
crystal. This may well be the cause of various 
hypersensitivity effects produced by mercury 
vapor. 

Similar reasoning leads to the conclusion that 
the replacement of some of the Br~ ions by S-~ 
ions, either within the crystal or on the surface, 
should also make the transition from B/ to B3 
energetically easier. For neutrality, there should 
be one extra positive ion (Ag+) for each Br~ ion 
so replaced, each extra ion helping to make B3 
more stable than BJ. 

One would, in fact, expect a sulfide ion (in the 
crystal interior) to tend to produce a B3 type 
structure with two (or perhaps more) extra Agt 
ions around it, in addition to those normally 
around each negative ion. This region would 
attract extra electrons (e.g., photoelectrons) to 





nearby positive potential holes. Perhaps the 
transition to this modified B3 structure would 
occur only if and when such extra electrons are 
available. The structure at and immediately sur- 
rounding the S~~ ion, according to this picture, 
can properly be described as a ‘“‘silver sulfide 
speck.”” (The AgeS and AgBr structures are not 
separated by a distinct boundary, but are co- 
herent; some of the silver ions are common to 
both.) The hypothesis under discussion thus 
furnishes a reasonable explanation for the marked 
sensitizing effect, discovered by Sheppard,° of 
sulfide compounds which are capable of producing 
silver sulfide during the making of the photo- 
graphic emulsion. 

The adsorption of polar dyes onto the surfaces 
of a silver bromide crystal has several effects 
which are pertinent to the present discussion: 
(1) By occupying surface sites which would 
otherwise be occupied by Br~ or other negative 
ions, each adsorbed dye molecule or ion produces 
a region of excess positive charge—as compared 
with other surface regions, at least. This is 
especially true if the adsorbed dye is a positively 
charged ion. (2) Because of the positive charge in 
this region, extra electrons (photoelectrons and 
developer electrons) are attracted to it. (3) The 
presence of extra electrons tends to transform 
this region to a B3 type structure. The composi- 
tion of the dye may also be such as to favor this 


‘ 


°S, E. Sheppard, Phot. J. 65, 380 (1925). 
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BI 


jan 2) ~ 
J ical yf 
+ 
Fic. 2. Projection, onto a {110} plane, of a crystal, the upper portion 


of which has the B3 structure, the lower portion the B/ structure. 
{111} planes are horizontal, in both structures. 


transformation. (4) Absorption of light by the 
dye molecule or ion releases an electron, either 
from the dye or from an adjacent Br ion, in the 
very place where it would be most effective in 
producing a Bl-to-B3 transformation. (5) In 
the vicinity of the adsorbed dye molecule or ion, 
the potential barrier hindering the approach and 
adsorption of developer ions is less than elsewhere 
on the crystal surface, because the net charge 
density is less negative than elsewhere.19 1! 

For several reasons, therefore, one would ex- 
pect adsorbed dyes to have a photosensitizing 
effect, as they do. Obviously, different dyes 
should differ as regards their effectiveness. 


RELATIVE DIMENSIONS 


Unit cells of the BJ and B3 structures contain 
' the same number (4 each) of Ag+ and Br- ions. 
The unit distance a) computed for B3 is about 
6 percent larger than the experimental value‘ of 
ao for B1. The calculated B3 value might be in 
error, perhaps, by 2 or 3 percent, but probably 
not more. 

The replacement of some of the Br~ ions by I- 
should tend to make the Bi and B3 dimensions 
a little more nearly equal. The presence of extra 
Ag* ions might have a slight effect in the same 
direction. Extra electrons in the structure, how- 


10R, M. Evans and W. T. Hanson, Jr. Phot. J. 77, 497 
(1937). 
11 T, H. James and G. Kornfeld, Chem. Rev. 30, 1 (1942). 








ever (especially with extra Ag* ions also present), 
would have a much more marked effect. For a 
definite, but unknown, concentration of extra 
electrons (or of extra electrons plus extra Ag* 
ions), a transformation from B/ to B3 could take 
place without any change of dimensions. It could 
occur even in the interior of a single crystal, 
without breaking up the crystal structure. 

At the surface of a crystal there is necessarily 
considerable strain, because the environment of 
the surface atoms is not the same as that of 
atoms in the interior. This surface strain would 
be relieved somewhat if the surface atoms were 
in a structural arrangement which, under similar 
circumstances, would be slightly less dense. 
Provided a suitable transition is possible between 
the two structure types (see below), we might 
therefore expect a B3 type surface on a Bl 
crystal. Other factors, especially the nature of 
the environment outside of the crystal (adsorbed 
ions, etc.), are also important, however. It may 
be that the surface of a crystal in a photographic 
emulsion contains regions of both types, BJ and 
B3, or that the surface is all of the B1 type but in 
certain areas (sensitivity specks) can easily 
change to the B3 type on the addition of extra 
electrons. 


THE TRANSITION FROM BI TO B3 


The distribution of bromide ions, considered 
by themselves, is precisely the same for a B3 and 
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for a B/ structure. If the dimensions are also the 
same, we may picture a single crystal as being in 
part Bl and in part B3, the Br- ion distribution 
being continuous throughout the whole, only the 
distribution of Ag* ions and extra electrons being 
different in the two parts. The structure in a 
portion of such a double-type crystal is repre- 
sented in Fig. 2. There is no sharp break between 
the two structures, one sheet of Br~ ions being 
common to both. Where the two structures come 
together, as well as within each structure, the 
arrangement of ions is energetically reasonable. 

For simplicity and clarity, the boundary be- 
tween the two structures in this figure follows a 
definite {111} plane. This need not actually be 
so; an irregular boundary would be quite possible. 
The boundary between the two structures could 
shift very readily, one atom at a time. An Ag* ion 
on the B1 side of the boundary need shift only 
slightly, as shown in Fig. 3, in order to take its 
place as a portion of the B3 part of the crystal, 
the shift involving only a very small activation 
energy (less than 10-" erg). Moreover, the 
presence of extra electrons in the B3 region would 
result in an attraction for the Agt ions on the B/ 
side of the border and would make this energy 
hump still smaller. 

Extra electrons in any region of a crystal (over 
the number required to balance any excess 
positive charge which may previously have 
existed there) would attract positive ions, causing 
a (slow) migration in that direction of any 
interstitial Ag+ ions present in the vicinity. From 
the calculations of the preceding paper, it follows 
that equilibrium concentration of interstitial ions 
in a B3 region is considerably greater than that in 
a Bl region, also that the rate of migration of 
such ions through a B3 structure is considerably 
faster than that through a B/ structure. At any 
rate, an extra electron would before long be 
neutralized by an extra Agt ion near it, occupying 
one of the close sites of negative potential. This 
process prevents the accumulation in a B3 latent 
image of an excess negative charge—which would 
repel additional photoelectrons. 

Even after ‘‘neutralization” by an extra Agt 
ion, an interstitial electron in a B3 structure 
would still attract (over small distances) other 
interstitial Ag+ ions and, in time, would be ex- 
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Fic. 3. Illustrating the nature of the shift made by a 
silver ion when it leaves the B/ structure on one side of a 
B1-B3 boundary to become part of the B3 structure on the 
other side. 


pected to hold a second such ion in the immediate 
vicinity. This region would then have a net 
positive charge, making it (like the original 
sensitivity speck) an attractive center for more 
photoelectrons or developer ions. 

Reciprocity law failure at high intensities may 
be attributed to the accumulation of photoelec- 
trons in a B3 region faster than they can be 
neutralized by Agt ion migration or faster than 
the potential resulting from the charge can be 
attenuated by growth of the B3 region at the 
expense of Bl. Reciprocity law failure at low 
intensities may be ascribed to slow diffusion of 
captured electrons from B3 regions, or to re- 
gression of the B3 region after neutralization by 
Agt ions of the charges of the excess electrons 
therein (or, even more, after the acquisition of 
more Agt ions than enough to neutralize the 
photoelectrons). 

The Herschel effect may be attributed to the 
photoelectric ejection of extra electrons in a B3 
latent image region. This would cause the B/—B3 
boundary to shift in such a direction as to de- 
crease the size of the B3 region or to eliminate it 
entirely. 


DEVELOPMENT 


The process of development of a silver halide 
grain consists essentially of the addition of elec- 
trons and the removal of halide ions. Both optical 
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Fic. 4. Illustrating the formation of a metallic silver crystal 
(A1 structure) on a {100} face of a B3 AgBr crystal. 


microscope studies! !* and electron microscope 
studies'*-!* have shown that in some cases the 
metallic silver produced by development grows 
outward from the silver halide crystal surface. 
The places at which these growths start from the 
crystal may reasonably be assumed to be latent- 
image regions. 

A B3 latent image region at the crystal surface 
can be assumed to have a lower potential barrier, 
resisting adsorption of electronegative developer 
ions, than a B/ surface region (covered by nega- 
tive ions) for several reasons: (1) There is 
probably little or no excess of negative ions on 
the B3 surface, as already noted. (2) At the 
surface, as in the interior, a B3 structure has sites 
of positive potential—other than those at the 
normal ion positions—which would attract close 
negative ions; this is not true of a BJ structure. 
(3) As a result of the tendency of each extra 
electron in a B3 structure to accumulate one or 
more extra Agt ions in the immediate neighbor- 


2 W. Scheffer, Brit. J. Phot. 54, 116, 271 (1907). 
13 C. Tuttle and A. P. H. Trivelli, Phot. J. 68, 465 (1928). 
4M. von Ardenne, Zeits. f. angew. Phot. 2, 14 (1940). 
ps C. E. Hall and A. L. Schoen, J. Opt. Soc. Am. 31, 281 
(1941). 
16F, F. Lucas, see reference 1, p. 834. 
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hood, each B3 region would (sometime after ex- 
posure) either be neutral or (more probably) 
positively charged. (4) The extra electrons—and, 
to a lesser extent, the extra Ag+ ions—in a B3 
region are very mobile, shifting readily on the 
approach of another charged particle. In other 
words, a B3 latent image is readily polarized on 
the approach of a developer ion, this polarization 
producing an attraction. A BI region is not 
readily polarized. 

On adsorption of a developer ion to a B3 sur- 
face, the former gives up an electron to the latter. 
Extra electrons entering the B3 lattice in this 
way attract silver ions from below the surface, 
pulling them over energy humps which are small, 
to previously unoccupied sites of negative po- 
tential in and near the surface. This produces a 
distribution of silver kernels in the {100} planes 
exactly like the distribution in a {100} plane of 
metallic silver, even as regards dimensions."" (See 
Fig. 4.) A surface of this sort may quite properly 
be described as a surface of metallic silver. 

Additional developer ions might be adsorbed 
to such a silver surface even more readily than to 
a B3 surface. Whether or not that is the case, 
continued action of the developer furnishes more 
electrons, which move toward this region of 
excess silver kernels. As the concentration of 
electrons builds up, a second sheet of silver 
kernels is formed, outside of the first, the ions 
going to the energetically most favorable sites, at 


17 The simple relationship existing between the dimen- 
sions of metallic silver and those of ordinary (B/) silver 
bromide has been noted by J. J. Trillat and H. Motz, 
J. de phys. et rad. [7] 7, 89 (1936), by P. D. Dankov, 
Comptes rendus U.R.S.S. 24, 733 (1939), by P. D. Dankov 
and A. A.'Kochetkov, ibid., 26, 785 (1940), and by H. 
Wilman, Proc. Phys. Soc. London 52, 323 (1940). Dankov’s 
assumed mechanism for the formation of a metallic silver 
crystal on the surface of a BIJ AgBr crystal, however, is 
more complicated and less reasonable than the mechanism 
postulated here for the production of a silver crystal on a 
B3 AgBr crystal. Wilman showed by electron diffraction 
methods that reaction of halogen vapors on thin silver 
films produced silver halide crystals which were oriented, 
relative to the underlying silver structure, in the manner 
indicated in Fig. 4. For simplicity of description it is 
assumed, in Fig. 4 and in the explanation (in the text) of 
the formation of a silver crystal on a silver halide surface, 
that the former starts at a {100} surface of the latter. That 
is not theoretically necessary, however. Oriented growth of 
a silver crystal on a {111} surface of a silver bromide 
crystal would seem to be especially likely, since the 
arrangement of silver kernels in a {111} face of a B3 (or B1) 
silver bromide crystal is just the same, even as rege urds 
dimensions, as that in a {111} face of a metallic silver 
crystal. 
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points opposite the centers of the squares of the 
underlying silver layer. They may reach these 
sites by migration along the crystal surface from 
points where they are least tightly held, or they 
may come from the previously formed silver 
layer, the sites vacated being soon refilled by 
other like ions migrating (because of the po- 
tential gradient) from the crystal interior. Under 
certain circumstances, silver ions might also be 
deposited from the solutior.'®"! 

Continuation of this process produces a silver 
crystal growing outward (possibly also spreading 
sideways, under suitable conditions) from the 
original surface, the two structures fitting to- 
gether as indicated in Fig. 4. Without sideways 
growth, a filamentary structure would thus be 
produced. 

The growth of the silver crystal is necessarily 
accompanied by diffusion of electrons, either 
from it or directly from the developer ions, into 
the B3 crystal, also by growth of the B3 region at 
the expense of the B/ part of the crystal. 

In view of the considerably higher activation 
energy for motion of a silver ion through a Bl 
structure than through a B3 structure,’ it seems 
probable that such motion occurs only in the B3 
regions or along a BJ (or B3) surface. The 
migration of the Agt ions continually diminishes 
the forces holding the surface bromide ions. 
These, therefore, become free to diffuse, more or 
less readily depending on the composition of the 
solution and the nature of the surface, into the 
solution. 

Very long exposure to light, in spite of the 
growth of the B3 regions and the migration of 
Agt ions, must eventually cause the negative 
charge on these regions to build up, eventually 
resulting in the formation and growth of metallic 
silver crystals, as when the extra electrons are 
supplied by the developing agent. Uncharged 
bromine atoms and, from them, bromine mole- 
cules are of course produced at the same time. 


SUMMARY 


The processes involved in latent image forma- 
tion and growth and in subsequent development 
of the crystal grain may be outlined as follows: 

1. Absorption of a light quantum by an atom 
(usually Br-) produces a photoelectron. 


PHOTOGRAPHIC 
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2. Unless recaptured by the atom from which 
it came (or some other uncharged bromine atom), 
the photoelectron diffuses through the crystal 
until it settles in a site having a positive potential 
somewhat greater than that at any point within 
the normal structure. Such sites (sensitivity 
specks) will usually be at the crystal surface, near 
ions or molecules of particular kinds adsorbed on 
its surface. Regions containing sulfide ions sur- 
rounded by extra silver ions (‘‘silver sulfide’), 


‘mercury ions with or without extra electrons 


(‘‘metallic mercury”’), extra silver ions with or 
without extra electrons (‘‘metallic silver’), and 
crystal surface areas covered by appropriate 
adsorbed dye ions (instead of the extra halide 
ions elsewhere on the surface) would be expected 
to be especially effective in attracting and 
trapping photoelectrons. 

3. Cooperation between the electron and these 
foreign ions or molecules produces a shift, over a 
limited region in the immediate vicinity, from the 
normal BJ type of structure to a B3 (or similar) 
type. The transformation involves only slight 
movement of the silver ions. 

4. Any net negative charge produced in a B3 
region by the addition to it of (one or more) 
photoelectrons is slowly neutralized by the migra- 
tion of silver ions toward these electrons. Further 
migration and accumulation of Ag+ ions around 
each electron may even give the region a net 
negative charge. Conditions are then favorable 
for the further addition of more photoelectrons or 
for development. 

5. Additional photoelectrons finding their way 
to this B3 region produce further shifts of silver 
ions at the B3—B1 boundary, thereby causing 
this boundary to shift in such a way as to increase 
the size of the B3 region. 

6. In the development process, extra electrons 
from the developer enter the surface of a B3 
region (more readily than that of a BJ region). 
After development has started, they may also 
enter at the silver surface. These electrons diffuse 
into the interior of the B3 region, causing it to 
expand at the expense of the B/ part of the 
crystal, as outlined in (5). 

7. The extra electrons also produce migration 
especially in the B3 structure and 
both 





of silver ions 


along the crystal surface, where the 
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underlying structure is B3 and where it is Bl— 
in such a way as to start a crystal of metallic 
silver on the surface. 

8. Further additions of migrating silver ions 
cause this silver crystal to grow. Because of the 
difficulty of concurrent migration of bromide ions 
away from the silver crystal, the latter usually 
grows outward from the original surface of the 
silver halide crystal, producing the filamentary 
structure observed. 

9. The ionic migrations just described leave 
the bromide ions on the surfaces held only very 
loosely or not at all; the excess negative charge 
within the crystal in fact repels them. They there- 
fore diffuse away into the solution, at a rate de- 
pending on the nature of the solution, on the rate 
at which electrons are being supplied to the 
crystal aggregate by developer ions, and on 


JR., 





AND J. W. McBAIN 
the rate of migration of silver ions within the 
aggregate. 

It should be emphasized that in none of the 
assumed processes is any break-up and _ re- 
formation of crystals necessary : both the transi- 
tion from B/ to B3 and that from B3 to metallic 
Ag are continuous. The B3 structure “catalyzes” 
the development process by making it possible 
for the addition of electrons from the developer, 
the formation and growth of silver crystals, and 
the removal of Br~ ions from the silver halide 
crystal all to occur at different, not necessarily 
adjacent, locations. 

In conclusion, the writer takes pleasure in 
acknowledging the help received from discussions 
of this subject with other members of the Kodak 
Research Laboratories staff, especially Drs. S. E. 
Sheppard, T. H. James, and J. H. Webb. 





THE JOURNAL OF CHEMICAL PHYSICS 


VOLUME 11, 


NUMBER 9 SEPTEMBER, 1943 


X-Ray Diffraction Investigation of Sodium Stearate from Room Temperature to the 
Melting Point 
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(Received March 11, 1943) 


Sodium stearate of the gamma-polymorphic form was photographed by Cu Ka x-rays in a 
powder camera at seven temperatures between room temperature and 301°C. The long spacing 
is constant up to 132°C, the transition temperature between the subwaxy and waxy phases, 
whereupon it continues to decrease markedly in length up to the melting point. The two most 
intense short spacings increase in length up to about the transition between the stable waxy 
crystalline phase to superwaxy phase at 167°C. Both lines then appear to coalesce to form a 
single line up to above the final melting point, from liquid crystal to liquid, at 288°C. The 
space lattice was assumed to be monoclinic over the whole temperature range studied, with 
hexagonal close packing of the chains themselves above 167°C, possibly as a result of molecular 


rotation of the chains about their axes. 


INTRODUCTION 


ODIUM stearate exhibits seven stable phases 

between room temperature and the melting 

point at 288°C.!:* The phase changes are given in 
Table I. 


The successive changes of sodium stearate on 
warming have been observed by the dilatometer 


* Present address: 64 Terrace Road, Medford, Massa- 
chusetts. 

1R. D. Vold, J. Am. Chem. Soc. 63, 2920 (1941). 

2M. J. Vold, M. Macomber, and R. D. Vold, J. Am. 
Chem. Soc. 63, 168-175 (1941). 





method,* the hot wire method of I. Fankuchen 
with crossed Polaroids, or with the ‘mikro- 
polychromar”’ microscope,? and the calorimetric 
method,! and it was thought the x-ray diffraction 
method would give additional information. 

The gamma‘~* polymorphic form was used in 

?R. D. Vold and M. J. Vold, 61, 808-816 (1939). 

4A. de Bretteville, Jr., and J. W. McBain, Science 96, 
470-471 (1942). 

5M. J. Buerger, L. B. Smith, A. de Bretteville, Jr., and 
F. V. Ryer, Proc. Nat. Acad. Sci. 28, 526-529 (1942). 


6 J. W. McBain, A. de Bretteville, Jr., and S. Ross, 
J. Chem. Phys. 11, 179-183 (1943). 
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this investigation. It may be obtained by heating 
the hydrated alpha-form in an open tube, where 
it changes to the beta-form at 54°C and then to 
the supercurd phase (stable between 89°C and 
114°C); thereafter the stearate reverts to the 
gamma-form, which is the most stable form at 
room temperature. The gamma-form exhibits 
at room temperature the same long x-ray spacing, 
even when it has been obtained by cooling 
sodium stearate from above the melting point. 
W. W. Lee has shown that sodium stearate forms 
a hemihydrate when exposed to half-saturated 
water vapor.’ However, all fibrous forms of soap 
sorb water in a continuous manner, depending 
upon the relative humidity. Hence any actual 
analytical composition (other than the hemi- 
hydrate) is devoid of significance. McBain, 
Boldnan, and Ross* have presented an x-ray 
examination of the hydrated gamma-form of 
sodium stearate. 


EXPERIMENTAL METHOD 


A circular x-ray heating camera’® of the pinhole 
type, was used to take pictures at successively 
increasing temperatures of the sample. The 
sodium stearate used was made from Kahlbaum 
acid by Dr. E. L. McBain.'!° The sample to film 


TABLE I. 











Phase Temperature Method 





curd-supercurd 90°C dilatometer, x-rays 
subwaxy 117°C dilatometer, x-rays 

waxy 132°C dilatometer, hot wire 
superwaxy 167°C ___ dilatometer, hot wire 
subneat 205°C ___ dilatometer, hot wire 

neat 257°C ___ dilatometer, hot wire 
isotropic [liquid ] 288°C dilatometer, hot wire, visual 

observation, microscope 








radius was 5.73 cm. The sample was sealed in 
clear radio code No. 707 Corning glass not over 
l-mm inside diameter and rotated at 1 r.p.m. 
Cu Ka radiation was obtained from a copper 
x-ray tube with a thin Pyrex “bubble window” 


7™W. W. Lee, Ph.D. Thesis, Stanford University, (1941). 
Obtainable in bibliofilm from the U. S. Department of 
Agriculture. J. W. McBain and W. W. Lee, Ind. Eng. 
Chem. 35, 784 (1943). 


’J. W. McBain, O. E. Boldnan, and S. Ross, J. Am. ° 


“hem. Soc. (in press). 

* Rev. Sci. Inst. 13, 481-483 (1942); J. App. Phys. 14, 
137 (1943), 

E. L. McBain, J. Phys. Chem. 44, 1014 (1940). 


STEARATE 


0 


Fic. 1. Short spacings in angstrom units as a function 
of temperature. The phases are bounded by the dashed 
lines, and are designated by: a(curd); 6(supercurd) ; 
c(subwaxy) ; d(waxy) ; e(superwaxy) ; f(subneat) ; g(neat) ; 
h(isotropic liquid). 


made by James F. Lee,** under a grant furnished 
through the Cyrus M. Warren Fund. The radia- 
tion was filtered by 0.01-cm foil of electrolytic 
nickel. The voltage applied to the tube was 28 kv 
at 28 ma. The exposure time was 2 hours. 

The photographs were measured visually with 
a sliding cross hair and vernier reading to 75°. 


EXPERIMENTAL RESULTS 


Figure 1 shows the short spacings d, as a 
function of the temperature. These data were 
obtained from the photographs shown in Fig. 2. 
Starting from room temperature the two most 
intense side spacings increase in length and in 
diffuseness as the temperature increases. This 
diffuseness may be due to decreased alignment 
because of the release of constraint with increasing 
temperature. Finally these two lines appear to 
merge above 170°C in the superwaxy phase 
(167°C-205°C). 

The long spacing d is tabulated in Table II as 
a function of temperature. The long spacing 
begins to decrease at 145°C in the waxy phase 
(132°C-167°C) and continues to decrease as the 
temperature increases. The value at 301°C in the 
isotropic phase is the mean of results of measure- 
ments on oriented doublets of the second and 
third orders. The doublets are short sharp arcs, 
in the 9 and 3 o'clock positions on the photo- 
graph, parallel to the axis of the sample capillary. 


** 193 Willow Road, Palo Alto, California. 








428 A. DE SBRETTEVILGSE., 


The oriented doubling may be due to the fact 
that a thin layer of molecules near the capillary 
walls is oriented. R. D. and M. J. Vold noticed 
that the soap becomes oriented in the isotropic 
liquid and neat phases between glass cover plates 
used in the hot wire method." 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 





Fic. 2. Successive changes of sodium stearate diffraction 

attern (@ range 0°-44°) with increase of temperature. 

he photographs show temperatures at: (a) room; (b) 
125°C in the subwaxy phase; (c) 145°C in the waxy phase; 
(d) 170°C in the superwaxy phase; (e) 220°C in the sub- 
neat phase; (f) 269°C in the neat phase; (g) 301°C in the 
isotropic liquid. 


" See reference 3, p. 171. 





L., BED J. W. McBain 
TABLE II 
ie d(A) n(order) Intensity 
25°C 44.3 5 s 
44.9 7 Ww 
4.68 m 
4.05 Ss 
2.95 Ww 
125°C 44.7 3 m 
4.83 WwW 
4.30 s 
2.97 Ww 
145°C 42.9 3 s 
4.81 m 
4.36 s 
2.99 WwW 
170°C 39.9 3 m 
4.68 m 
2.98 vw 
220°C 37.9 2 s 
38.3 3 Ss 
4.78 m 
269°C 34.0 2 s 
33.3 3 m 
4.92 m 
301°C 31.9 2 s 
32.3 3 m 
4.82 WwW 
DISCUSSION 


Powder diffraction data are given in references 
4 and 6, and the fiber diagram with derived 
unit cell in reference 8 for the gamma-form of 
sodium stearate at room temperature. The alpha- 
lattice was investigated by M. J. Buerger with 
single hydrated crystals of sodium stearate and 
sodium palmitate grown by A. de Bretteville, Jr., 
and F. V. Ryer.” It has monoclinic symmetry. 
Crystals of the beta-form obtained by heating 
the alpha to 54°C are shattered and opaque and 
poorly suitable for analysis. P. A. Thiessen and 
J. Stauff characterized the beta-form as mono- 
clinic,"* and McBain, de Bretteville, and Ross 
used their rotation photograph to index correctly 
the beta-powder photograph. 

As shown in Table II, the long spacing 
decreases as the temperature increases, especially 
above about 132°C. This may be attributed to 
an increase of tilt of the chains. The length of 
the ¢ axis is uncertain, because for a monoclinic 


2M. J. Buerger, Proc. Nat. Acad. Sci. 28, 530 (1942). 
13 Zeits. f. physik. Chemie A176, 397 (1936). 
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space lattice the long spacing is c sin 8, where c, 
the length of the ¢ axis, and 8, the angle of 
inclination of this axis with respect to the basal 
plane, are both unknown. However, if the double 
molecules of length 52.1A lie in the ¢ direction, 
the value of 8 would be about 60°. Then, taking 
also the increase in side spacing into account, 
the x-ray data would confirm the dilatometer 
results that as the temperature is raised the 
specific volume of each successive phase is larger. 

Hexagonal cylindrical packing inclined at an 
angle to the basal plane might result from 
spinning of the carbon chains about their axes 
at high temperatures in the phases above 167°C. 
The subwaxy and the waxy phase are definitely 
crystalline, but the superwaxy phase, existing 
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between 167° and 208°C, more nearly resembles 
the liquid crystalline forms subneat (208°-238°) 
and neat (238°-288°). It is probable that in 
these three ferms the dipole carboxylate groups 
are fixed with respect to each other while the 
paraffin chains rotate on their axes. The inclina- 
tion of the axes must result from the fixation of 
the carboxylate heads which largely determine 
the angle of the single bond holding the paraffin 
chain. 

It may be noted that the analogous normal- 
paraffins also exhibit numerous phase changes 
below their m.p. with hexagonal cylindrical 
packing above the m.p., as shown by A. Miiller.™ 


4 A. Miiller, Proc. Roy. Soc. A138, 514-530 (1932). 
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An Absolute Method for the Determination of 
the Area of a Fine Crystalline Powder 


WILLIAM D. HARKINS AND GEORGE JURA 
University of Chicago and Universal Oil Products Company, 
Chicago, Illinois 
July 12, 1943 


RIOR tothe development of the methods presented here 
in three letters, there has been only one reliable method 
for the determination of the area of a solid in the form of 
thin sheets, of a powder, or a porous material. That 
method, developed by Emmett, together with Brunauer 
and Teller,! is based on an isotherm (BET) which does not 
give an area, but only the volume (v,,) of gas adsorbed 
upon the surface of the solid to form a complete monolayer. 
From this, the number (NV) of molecules in the monolayer 
may be calculated. From the density of liquid nitrogen, 
they calculate the area per nitrogen molecule as 16.2A? 
and irom the solid 13.6A?2, but neither of these can be the 
correct value to use. One of our objects is to determine 
the correct values for the different types of solids. 

The area of a very finely divided solid in the form of a 
powder may be determined by a calorimetric method. 
Due to the fact that our work has been restricted almost 
entirely to those porous solids of the greatest practical 
importance at the present time, it has been possible to 
apply the absolute method to only two materials thus far. 
These are: 


(1) a fine crystalline TiOz in the form of anatase, and, 
(2) a sample of the same material coated with a small 
amount of Al.O3. 


That the surface of the latter material is different from 
that of the former is shown by the fact that the presence 
of the Al,O; increased very considerably the energy of 
immersion (£;) of the material (per unit area) in water. 

In order to determine the area of the powder, it is 
outgassed in a Pyrex tube, connected to another tube 
which contains a liquid. A smaller tube, which connects 
those which contain the solid and the liquid, has a trans- 
verse glass membrane, which prevents the vapor of the 
liquid from reaching the solid. 

The tube which contains the solid is now sealed off from 
the vacuum system and placed in a good thermostat 
(e.g. at 25°C). The membrane is broken and the vapor 
of the liquid is allowed to come into equilibrium with the 
surface of the solid. If a large volume of solid is used, 
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glass tubes should lead the vapor to every part of the 
powder. The attainment of equilibrium may require a 
period of days. 

The principal condition to be met is that the liquid 
chosen for the determination must give a zero angle of 
contact with the solid. For example, with an unoxidized 
carbon, a hydrocarbon may be used, while with polar 
solids water is suitable. 

The next step is to immerse the powder, whose particles 
are coated with a film of the liquid, in a considerable 
volume of the same liquid in an extremely sensitive calo- 
rimeter. The amount of heat (g) in ergs per gram, divided 
by the total surface energy of the liquid (118.5 erg cm~* 
for water) gives the area per g of powder. Actually there 
is a correction of the order of one percent for the thickness 
of the film of liquid. 

This method is based upon the idea that the surface of 
such a duplex film as is formed by adsorption from satu- 
rated vapor, and with a zero contact angle, has exactly 
the surface energy of the liquid itself. With water, the 
weight of the film indicated a thickness of about 9 mo- 
lecular layers, while our measurements of energy show 
that the effect of the solid upon a water film cannot be 
detected beyond five layers by a calorimetric measurement. 

In order to obtain much accuracy by this method, it is 
essential to have powders with high areas. The powder of 
highest area gave 13.8 sq. meters per g, while by the use 
of Emmett’s usual assumption of 16.2A? for nitrogen, the 
BET method also gave 13.8 sq. meters per g. 


1S. Brunauer, P. H. Emmett, and E. Teller, J. Am. Chem. Soc. 60, 
316 (1938). 





A New Adsorption Isotherm Which Is Valid 
Over a Very Wide Range of Pressure 


GEORGE JURA AND WILLIAM D. HARKINS 
University of Chicago and Universal Oil Products Company, 
Chicago, Illinois 
July 12, 1943 


IL films on water exist as monolayers in a number of 
states, such as a two-dimensional perfect and im- 
perfect gas, the liquid expanded state, and several con- 
densed states. The film adsorbed on solids, although not 
restricted to monolayers, should also exhibit such states. 
Thus we have found that, in general, the very dilute films 
are monolayers in the gaseous state, as exhibited by their 
pressure (7) —area (c) relations. 
Condensed monolayers on water are characterized by 
the fact that they exhibit the linear pressure-area relation 


x=b—ae (1) 


in which a and 6b are constants and o is the mean area 
per molecule. 

The pressure-area relations of films on solids had not 
been determined, but when we plotted them for crystalline 
powders and porous solids, it was seen that Eq. (1) is 
valid for condensed films on solids. 

The remarkable difference is, however, that the con- 
densed film is monomolecular for only low values of the 
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film pressure, but the linear relation persists up to higher 
pressures where the film is several molecules thick. 

The straight lines were found to be parallel, or almost 
parallel, so the value of a is constant, or nearly constant, 
while that of 6 varies. 

Usually there is only one straight line (only one value 
of a), but in the case of some powders and porous solids, 
there is a second straight line at higher pressures, which 
represents a film of greater compressibility. Ona charcoal, 
where the film is almost monomolecular, the compressibility 
became much smaller at the higher pressures. 

Considerable calculation is involved in obtaining the 
values of the film pressure, so Eq. (1) was transformed into 
the equivalent relation 


log p=B-—A/v? (2) 
or, 


log (p/po) = B—A/v*, (3) 


which involves only the quantities » and v which are 
measured directly in the experimental determination of 
adsorption. 





ce Np adsorbed per g 














p-mmHg 


Fic. 1. New and simple BET isotherms of Nz on TiO: (anatase). 


This equation is remarkable in that (1) it is the simplest 
adsorption isotherm thus far found; (2) it is valid over 
more than twice the pressure range of any other two 
constant adsorption isotherms. 

Thus Fig. 1 shows that with TiOz in the form of anatase, 
our new isotherm agrees with the experimental data to 
both lower and much higher values of pressure than the 
BET isotherm. At very high pressures, both isotherms 
give too high values for the amount of gas adsorbed, but 
our curve lies very much closer to the experimental curve. 

With porous solids, Fig. 2, our equation seems at low 
pressures to be somewhat the better of the two, but gives 














THE EDITOR 431 
T T T T T T 
300 / 4 
PY ° 
—— NEW / 
osteo BET J / 
f = 
{oe 
200 a 4 
ao 
- fr 
“ 
100 f 4 
b 
4 
1 1 l l l l 
% 500 


p-mm ig 


Fic. 2. New and BET isotherms of N2 on a porous solid. 


too high amounts, while the BET method gives too small 
amounts, adsorbed. At higher pressures, our isotherm is 
much the better, but both isotherms give too high values 
at sufficiently high pressures. 

The obvious explanation for this particular deviation of 
our isotherm is that the smaller capillaries fill up,’ and 
prevent the further adsorption which would occur if the 
pores were sufficiently large. However, the BET isotherm 
gives too high values at the higher pressures even if pores 
are absent. 





An Adsorption Method for the Determination 
of the Area of a Solid without the Assumption 
of a Molecular Area, and the Area Occupied by 
Nitrogen Molecules on the Surfaces of Solids 


WILLIAM D. HARKINS AND GEORGE JURA 
University of Chicago and Universal Oil Products Company, 
Chicago, Illinois 
July 12, 1943 


N the preceding letter,! it was shown that the adsorption 
isotherm 


log p/po=B—A/v* (1) 


is valid over an extremely large range of pressures. This 
equation represents an S-shaped isotherm. 

If log p/po is plotted on the y axis (Fig. 1), (semi- 
logarithmic paper), a straight line of slope —A is obtained. 
It is very simple to show that the square roots of the 
slopes A are proportional to the areas of the solids. Let A 
be represented by s. Then the area (2) for one gram of 
solid is given by: 


r=ks}=4.06s!, (2) 
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Fic. 1. Adsorption isotherms of nitrogen on some porous solids. 
These isotherms represent the surfaces of porous solids which have 
areas, beginning at the top, of 321, 365, 395, 409, 438, and 455 square 
meters per gram. These areas were determined from these slopes on a 
large scale plot. : 


in which & is the area in sq. meters per gram if v in Eq. (1) 
represents the volume adsorbed by 1 g of material. 

The value of the constant (4.06) was determined by the 
substitution in Eq. (2) of the value 13.8 sq. meters per 
grant for = as determined for a sample of anatase by the 
absolute method described in the first letter. The value of 
s was determined from the adsorption isotherm of this 
same powder. 

Although the isotherm of Brunauer, Emmett, and Teller 
is very different from the new isotherm, and although the 
two were developed from extremely different points of 
view, the agreement between the values for the area as 
obtained by the two methods is remarkable. For example, 
nitrogen was adsorbed on 56 porous and 4 non-porous 
solids, and the areas were calculated by the use of both 
isotherms. Since Emmett uses either 13.6A? or 16.2A? as 
the area per nitrogen molecule, let us choose a value 
intermediate between the two: 15.25A?. If this is done, 
it is found that for 58 of the sixty solids, the two methods 
agree to better than +9 percent. 

A more natural method of comparison is obtained if it 
is remembered that the BET method does not give area 
(=) but only the number (NV) of molecules to form a 
complete monolayer, while our method gives the area, but 
not the number of molecules. So, 


o=2/N (3) 


gives in each case the mean area (c) per nitrogen molecule. 
It is obvious that if this area o is used in the calculation 
of area by the BET method, the two methods must, of 
necessity, give identical results. 

Thus, from their isotherm, N may be obtained and 
from ours, 2, and from the two, the area per molecule of 
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TABLE I. Area occupied by nitrogen molecules on the surfaces of solids. 


Number of solids with No» areas 





between the limits Area per Ne 

1 13.6-13.7 
4 13.8-13.9 
6 Peak 14.0-14.1 
3 14.2—14.3 
2 14.4—14.5 
1 14.6—14.7 
3 14.8-14.9 
7 15.0-15.1 
9 Peak 15.2-15.3 
6 15.4-15.5 
2 15.6-15.7 
3 15.8-15.9 
4 16.0-16.1 
5 Peak 16.2-16.3 
3 16.4-16.5 
0 16.6—16.7 
1 16.8-16.9 
1 21.3 

61 


nitrogen on any solid may be calculated. The results of 
determinations for sixty-one solids are given in Table I. 

All of these solids were porous, except four, which gave 
areas for the nitrogen molecules as follows: anatase 16.1; 
anatase, sample two, 16.2; anatase coated with Al.O3, 16.2; 
zirconium silicate, 16.4. All of these gave nitrogen areas 
practically equal to that (16.2) calculated from the density 
of liquid nitrogen. Other determinations for crystalline 
powders should be obtained before it is concluded that 
this is true in general. 

It should be noted that only five values fall outside the 
range of areas calculated from solid (13.6) and liquid 
(16.2) nitrogen, and three of the five have values of 
16.4—16.5 which are almost within this range. It is possible 
that the single high value (21.3) is erroneous. 

It may be noted that the values in Table I give a 
distribution curve with three peaks, with the highest at 
an area of 15.25A, the second highest at 14.05, and the 
lowest at 16.25. 

If the relation between log p and 1/v? is expressed by 
two straight lines, then that for the Jower pressures is 
always taken, since this is the one in which the transition 
from a monolayer to a polylayer always occurs. 

Some porous charcoals have such fine pores that the 
nitrogen films are restricted to almost a single layer of 
molecules. For the only charcoal thus far used, our method 
agreed with that of Langmuir, which gave 1350, while 
the BET method gave 1050 sq. meters per g. In all other 
cases, our method agrees with the BET method, while 
Langmuir’s isotherm was found to be valueless, since it 
applies only to monomolecular adsorption. 


'G. Jura and W. D. Harkins, J. Chem. Phys. 11, 430 (1943). 





The Configuration of 1,3-Butadiene 


R. S. Rasmussen, D. D. TUNNICLIFF, AND R. ROBERT BRATTAIN 
Shell Development Company, Emeryville, California 
July 22, 1943 


ULLIKEN!? has calculated the intensity distribution 

in ultraviolet absorption (1700-2300A) for the cis- 

and trans- forms of 1,3-butadiene. Comparison with experi- 
mental data* indicated that at room temperature the 
trans- form predominates. Rieke and Mulliken? made 
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further calculations on the relative energies of the ground 
state of the cis- and trans- forms and concluded that 
trans-butadiene was the lower energy form, although 
Mulliken indicated that the validity of this conclusion is 
doubtful because of the approximations used. 

Two observations of interest in this connection have 
been made in these laboratories. 

(1) It was found that the ultraviolet absorption of the 
vapor in the region 2200-2400A is strongly temperature 
dependent, the optical density increasing with temperature 
(Fig. 1). The most probable explanation of this effect is 
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Fic. 1. Optical density of 1,3-butadiene, 1-cm cell. 


the existence of appreciable amounts of both isomers at 
room temperature and a rapidly changing equilibrium in 
this temperature region. Unfortunately, it was impossible 
to cover a wide enough temperature range to determine 
the three parameters governing this temperature depend- 
ence, namely, the absorption coefficients of each of the 
two forms and the energy of isomerization. 

(2) A redetermination was made of the _ infra-red 
absorption spectrum of the vapor using better resolution 
and more accurate wave-length calibration than had been 
used in previously published work.** The Raman spectrum 
of the liquid was obtained by Bradacs and Kahovec® and 
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was presumably taken at the temperature of a dry ice 
bath. Of ten infra-red bands and nine Raman lines observed 
between 1500 and 300 cm.-', six frequencies coincide 
within the sum of the experimental errors of the two 
determinations, as shown in Table I. Three of the coinci- 
TABLE I. Infra-red and Raman frequencies <1500 cm™ 
of 1,3-butadiene. 


Raman (liquid) 
(Reference 6) 


Infra-red (vapor) 
(present work) 


1488(1) 

1432(1)— —1437(4)P 

1383(5) 

1302(1)—— -—— ---— —1303(2)P 

1273(3)—— - ~— = —1276(4)P 
1200(4) P 

1014(9)— —1025(0) 

988(2) 

909(10)- —914(2) | 
899(1) | 

665(1)* 

520(5)*——__________———-513(4)P 
340(00) 


* Taken from reference 4. Raman polarizations and approximate 
Raman and infra-red intensities shown. 


dences occur with the strongest Raman lines in this range. 
Since the trans- form, due to its center of symmetry, 
should show no coincidences, it is indicated that, under 
the low temperature conditions of the Raman work, the 
cis- form is predominant. A breakdown of selection rules 
in the Raman case is not involved, since four coincidences 
occur with polarized, and hence allowed, Raman fre- 
quencies. 

From the foregoing observations and considerations it 
appears that cis-butadiene is the lower energy form. The 
cis- form predominates at the temperature of a dry ice 
bath, while at room temperature the trans- form predomi- 
nates, as indicated by the ultraviolet spectrum, although 
a considerable proportion must be in the cis- form in order 
to account for the coincidences of the infra-red spectrum 
observed at room temperature and the Raman spectrum 
observed at the lower temperature. 

1R.S. Mulliken, J. Chem. Phys. 7, 121 (1939). 

2R. S. Mulliken, Rev. Mod. Phys. 14, 269 (1942). 

3 E.g., see E. P. Carr, L. W. Pickett, and H. Stiicklen, Rev. Mod. 
Phys. 14, 260 (1942). 

4E. Bartholomé and J. Karweil, Zeits. f. physik. Chemie B35, 448 
eB. Barnes, U. Liddel, and V. Z. Williams, Ind. Eng. Chem. 


Anal, Ed. 15, 83 (1943). 
6 K. Bradacs and L. Kahovec, Zeits. f. physik. Chemie B48, 63 (1940). 
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